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ABSTRACT 


Modern phased array radar uses multifunction subarray antennas in a distributed fashion. 
Distributed subarrays (DSA) have the advantages of more efficient scheduling of track 
and search functions, rapid steering capability, decreased complexity in digital 
beamforming and better angular resolution. However, one disadvantage of the DSA are 
the extra grating lobes due to large subarray spacing which can cause ambiguities in 
angle measurements and excess background clutter. A possible approach to suppress the 
grating lobes is to design separate transmit and receive subarray antennas that have 
different radiation patterns. 

The purpose of this research was to develop a program based on the principle of 
pattern multiplication to synthesize and access the two-way antenna pattern for DSAs. 
The program, written in MATLAB, allows the user to study the two-way antenna pattern 
for different subarray architectures. The program was able to synthesize the pattern for 
isotropic elements, short dipoles and half-wave dipoles in a planar array above a ground 
plane. A simulation tool was also developed to map the grating lobe and null locations of 
the antenna patterns in direction cosine space. Several DSA configurations were 
examined, and the results showed that undesired grating lobes can be suppressed by 
subarray nulls. 
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EXECUTIVE SUMMARY 


Modern phased array radar uses multifunction subarray antennas in a distributed fashion. 
Distributed subarrays (DSA) have the advantages of more efficient scheduling of track 
and search functions, rapid steering capability, decreased complexity in digital 
beamforming and better angular resolution. One disadvantage of the DSA are the extra 
grating lobes due to large subarray spacing which can cause ambiguities in angle 
measurements and excess background clutter. A possible approach to suppress the grating 
lobes is to design separate transmit and receive subarray antennas that have different 
radiation patterns. 

The primary objective of this research was to develop a simulation tool to 
investigate the behavior and effectiveness in suppressing undesired grating lobes using 
the approach of two-way antenna pattern in DSA design. The fundamental array theory 
and principle of pattern multiplication, which fonned the basis of the simulation tool 
design, were discussed. With the principle of pattern multiplication, grating lobes can be 
suppressed by placement of subarray nulls at grating lobe locations. A simple program 
was developed in MATLAB to allow the user to visualize the placement of grating lobes 
and nulls in direction cosine space for a DSA configuration in the visible region. A 
simulation tool with graphical user interface (GUI) was developed and implemented in 
MATLAB to perform the two-way antenna pattern and power gain calculations for user 
configured DSAs. The program is capable of perfonning two-way pattern and power gain 
calculations for linear or planar DSAs consisting of isotropic elements, half-wave dipoles 
or short dipoles above a ground plane. The program is able to present the simulation 
results in the 9 pattern cut, (f> pattern cut, or as a three-dimensional mesh plot in 
direction cosine space. The program GUI provides a convenient way for the user to tweak 
the design configurations very quickly by changing the DSA parameters. 

The effectiveness of the two-way pattern multiplication approach to suppress 

undesired grating lobes by placement of subarray nulls at the grating lobe locations was 

demonstrated using the simulation tool. For a transmit DSA with subarray spacing of five 

wavelengths and unifonn amplitude illumination, a sidelobe level (SLL) of-25.5 dB was 

xvii 



achieved. Using the simulation tool, we demonstrated that low SLL and narrow half¬ 
power beamwidth (HPBW) of the two-way antenna pattern can be achieved using thinned 
transmit and receive arrays consisting of widely-spaced subarrays with non-coincident 
grating lobe locations. The simulation results showed that a peak SLL of -49.7 dB and 
HPBW of 0.9° were obtained for broadside illumination. Simulation was also carried out 
to examine the effects on the two-way pattern when different scan angles were applied. 
The simulation results showed that beam broadening was evident and peak SLL was 
increased for the two-way antenna pattern when large scan angles were applied to the 
main beams of the transmit and receive arrays. It was assumed that the main beams of the 
transmit and receive arrays were scanned to the same angle. 

The simulation tool that was successfully developed for this research will serve as 
a useful tool for both students and electromagnetic professionals to detennine and study 
the two-way pattern and power gain of different transmit/receive DSA designs. The 
simulation tool has also demonstrated the advantages of using separate transmit and 
receive antenna patterns to suppress undesired grating lobes and achieve narrow two-way 
beamwidth using fewer antenna elements. 



ACKNOWLEDGMENTS 


I would like to express my sincere gratitude to my thesis advisor, Professor David 
C. Jenn for his patience, support and advice. Without his valuable knowledge and 
guidance, this thesis would not have been completed. 

I would like to thank Professor Tri. T. Ha, for agreeing to be the second reader of 
this thesis. 

I would like to thank my lovely wife, Justina Chua, for her unwavering support, 
understanding and patience, since I could not spend enough time with her while I was 
preparing for this research. 

Lastly, I would like to thank Singapore Technologies (Electronics), and my 
superiors, Mr. Teai Yam Koon and Mr. Yeo Eng Choon, for providing me this 
opportunity to study for my master’s degree at the Naval Postgraduate School. 


xix 



THIS PAGE INTENTIONALLY LEFT BLANK 


xx 



I. INTRODUCTION 


A. BACKGROUND 


A phased array is a directive antenna made up of a number of individual antennas, 
or radiating elements [1]. These individual antennas are geometrically arranged and 
excited by relatively phased currents to produce desired radiation patterns. Since the 
1980s, phased arrays have gained wide use in military radar and communication 
applications. The Aegis ships in the US Navy use phased arrays for the AN/SPY-1 
multifunction radar (MFR) [2]. The phased array for the AN/SPY-1 radar is shown in 
Figure 1. 



Figure 1. AN/SPY-1 phased array radar (From [3]). 


With the advent of stealth technology to reduce the radar cross section (RCS) of 
modem naval fighting ships, platform real estate has become a precious commodity. This 
leaves limited areas on the ship structure for large antenna arrays. This is a challenge for 
design of antenna apertures for shipboard radar systems. One design approach is to use 
multi-function subarray antennas in a distributed fashion. A distributed subarray (DSA) 
consists of multiple subarrays of antenna elements physically distributed over a platform 
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to form a long baseline, very thinned array for accurate angular location of targets [4]. 
Synthetic apertures can be generated using coherent processing to produce optimized 
patterns for specific radar applications. DSAs have the advantages of more efficient 
scheduling of track and search functions, rapid steering capability, decreased complexity 
in digital beamforming, and better angular resolution. One disadvantage of the DSA is 
the extra grating lobes that are introduced due to large subarray spacing of multiple 
wavelengths. Grating lobes, a form of undesired aliasing, exist when more than one 
period of the array factor appears in the visible region. Grating lobes are undesirable 
because they have the same magnitude as the main beam, which can cause ambiguities in 
angle measurements and excess background clutter. 

A possible approach to suppress the grating lobes is to design separate transmit 
and receive subarray antennas. In radar systems, the received power depends on the 
product of transmit and receive antenna gains. In theory, the undesired grating lobes can 
be suppressed by combining transmit and receive antenna pairs that have different 
radiation patterns. The separate receive antenna can be designed with nulls in the 
direction of the transmit antenna array’s grating lobes in order to cancel or reduce its 
effect by taking advantage of the two-way antenna gain in the radar range equation. 

The major advantages of using separate transmit and receive antenna patterns are 
the following: 

• A narrow two-way beamwidth can be achieved using fewer antenna 
elements. 

• An increase in transmit/receive isolation can be achieved. 

• There is more flexibility in the physical placement of the antennas. 

• There is an added dimension of beamforming for digital phased arrays. 

• There is a possibility of clutter reduction using pattern design. 
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B. PREVIOUS RESEARCH 


The concept of distributed array radar (DAR) and factors involving two-way 
effective gain patterns for DAR arrays are discussed in [4]. The study described a 
convolution approach to evaluate DAR resolution patterns using different transmit and 
receive array apertures. Lin [5] described methods of improving angular resolution of 
shipboard radar using DSAs. His research also proposed using two-way patterns to 
suppress grating lobes in widely spaced subarrays. 

C. SCOPE OF RESEARCH 

The objective of this research is to design a simulation tool to investigate the 
behavior and effectiveness in suppressing undesired grating lobes by using the approach 
of two-way antenna pattern in DSA design. The element factor for short dipoles and half¬ 
wave dipoles above a ground plane is included. The array factor and principle of pattern 
multiplication are employed to reduce the computational burden. A program is developed 
using MATLAB to synthesize and evaluate the two-way antenna pattern for planar 
phased arrays. Several DSA configurations are examined, and the synthesis results 
presented and discussed. 

D. ORGANIZATION OF THESIS 

The overall thesis report consists of five chapters. The scope of the thesis and 
background, which includes a brief introduction of distributed phased arrays, is covered 
in Chapter I. 

The fundamental theory of phased arrays is discussed in Chapter II. The scanning 
pattern for arrays with elements above a ground plane is presented in this chapter. The 
array factor, ground plane factor, and element factor are introduced. In addition, the 
effects of grating lobes due to periodic variations in the array are explained. 

The principles of digital array beamforming, directivity, gain, and radar range 
equations for DSAs are introduced in Chapter III. The approach to suppress undesired 
grating lobes using the principle of pattern multiplication is also discussed. 
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The program developed and implemented in MATLAB for simulating the two- 
way antenna pattern for DSAs is described in Chapter IV. The simulation results 
conducted for several DSA configurations are presented and analyzed, as well. 

Finally, the summary and recommendations for future research and study are 
provided in Chapter V. 

The MATLAB codes for the program developed for this research are provided in 
the Appendix. 
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II. ARRAY PATTERN THEORY 


In this chapter, the fundamental theory of array antennas for a two-dimensional 
array above a ground plane is presented. The array pattern, array factor, element factor, 
and ground plane factor are introduced. Grating lobes in DSA design and their effects are 
also described. 

A. ARRAY ANTENNAS 

An array can be considered as a collection of identical radiating elements that are 
excited to achieve some desired radiating pattern. A linear array consists of antenna 
elements arranged in a single dimension (straight line). A linear array of dipoles is shown 
in Figure 2. A standard spherical polar coordinate system is used, with the elements along 
the x axis. Broadside is 0 = 0° and endfire is 0 = ±90°. 

z To far-field observation point at (r, 6,$ = 0) 



B. TWO-DIMENSIONAL ARRAYS WITH GROUND PLANES 

A two-dimensional (2D) or planar array consists of antenna elements arranged on 
a two-dimensional plane. The theories and formulas for the development of the 
MATLAB program to synthesize the two-way antenna pattern for planar DSAs are 
discussed in this section. The material can be found in [7]. 
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1. Array Factor 

The array factor (AF) is a function of the array geometry and relative excitations 
of the array elements. Consider a planar array of point sources laid out on a rectangular 
lattice in the x-y plane as shown in Figure 3. 



Figure 3. Planar array centered at the origin (From [7]). 


The pattern angles 6 and <j> for the planar array oriented in the x-y plane are 
defined by the coordinate system shown in Figure 4. 


2 A 



Figure 4. Array orientation and pattern angles 0 and (j) (After [8]). 

6 















If the array is centered at the origin, then the element locations can be written as 


and 


x„ = 


l„-(N x+ l) 


d X 1 * * * ? N v 


2m-(N y + \) 


y m ^ dy, tti !>•••> N y 


( 1 ) 


( 2 ) 


where d x and d y are the element spacings and N x and N y are the number of elements 
along the x and y axes, respectively. The array factor can be expressed as 


TV, N y 

AF= SI 

n= 1 m =1 


4 JPx„ sin 9 cos <j> JPy m sin 6>cos </> 
^mn c e 


( 3 ) 


where /? = 2n / A , and A is the wavelength. 

If the feeding arrangement results in a separable distribution (which is common), 
then A mn — A in A n , and 


AT 


AT, 


AF = AF x AF y = X A n e jpX " sin0C0 ^ X A m e j ^ m sin0C0 ^. 

n —1 m= 1 


( 4 ) 


Furthermore, if phases are introduced to scan the beam to the direction (6* s , </> s ) , then 


AT, 

I 

n =1 


N,. 


AF = x a e jfi x n{ sin0cos 0- sin0 s cos( / , s) X A ( sin ^ cos ^ _sincos A ) (5) 


m=l 


Finally, if the array is uniformly excited (A mn =1), the sum becomes a 
geometric series with the following closed form result [9]: 


AF = AF x AF y = 


sin(N x¥ J2)sm{N yVy l2) 


sm(y/ x /2) sin(i// y /2 


( 6 ) 


where 
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y/ x = J3d x (sin 6 cos (j> - sin 0 S cos </) s ) = J3d x (u - u s ) 
y/ y = fidy (sin 6 sin (f> - sin 0 S sin </> s ) = J3d y (v - v s ) 


( 7 ) 


and 

u = sin^cos^ 
u s , = sin 6 S cos <f> s 

( 8 ) 

v = sin^sin^ 
v v = sin 0 S sin (j) s . 

2 . Ground Plane Factor 

Ground planes are used typically to increase the directivity of the antenna array. 
Consider the array to be placed at a height h above a ground plane. If the ground plane is 
a perfect conductor of infinite extent, the method of images can be used to compute the 
field above the ground plane. The total array, including both the sources and their images, 
can be viewed as a planar array, where each element is a linear array along the z-axis as 
shown in Figure 5. 



Figure 5. Two-element linear array along the z-axis, where each element is a planar 

array (From [7]). 

With the images out of phase with the sources, the ground plane factor (GF) for 
the two element linear array along the z-axis can be written as 
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GF = e 


- o jPhcos0 _ -jph cos 0 _ 


= 2 j sin (/?/? cost?). 


(9) 


3. Element Factor 

The radiation behavior of a single antenna element is signified by its element 
pattern and can be characterized by its element factor. In general, the element factor has 
both 0 and (f> components. The element factors can be derived from the direction cosines 
(see Appendix A). 

For a half-wave dipole along the v-axis with a maximum feed point current I m , 
the element factors are 


and 


EFn = 


Pic, Pn e 


~ i P r 


2 nr 


cos 


V 2 


sin 0 cos (j) 


1-sin 2 t?cos 2 <f> 


cos 6 cos (j) 


EFa = - 


PlcJ nP 


~jP r 


2 nr 


cos 




— sin t? cos ^ 

v2 


1-sin 2 i?cos 2 (j) 


sin^ 


( 10 ) 


( 11 ) 


where rj 0 = 377 Q is the intrinsic impedance of free space. For a half-wave dipole along 
the y-axis, the element factors are 


and 


EFn = 


jVcJnP 


~jpr 


2 nr 


cos — sin 0 sin ^ 

V 2 _y 

1-sin 2 9 sin 2 (j) 


cos 6 sin (/) 


EF^ = 


iVo hn e 


~jpr 


2 nr 


cos 


sin 6 sin (j) 

J 'cos^. 


1-sin 2 i?sin 2 (j) 


( 12 ) 


( 13 ) 


For a half-wave dipole along the z-axis, the element factors are 
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(14) 


and 


EF n = 


j T Jo I m e ]Pr 

cos 

(Tcos*) 

2 nr 

sin 6 


EF+= 0 . 


(15) 


To normalize the element factors, the leading factor j?j 0 I m e nr is simply 

removed. If the dipoles are ideal (Hertzian with constant current and length of the dipoles 
L« A) instead of half-wave dipoles, then the terms in square brackets [ ] become 1, and 
L is added to the leading factor. Therefore, after normalization of the element factors, 
only the trigonometric functions to the right of the bracketed terms remain. 


4. Array Pattern 

The total nonnalized pattern factor of the array is obtained using the principle of 
pattern multiplication: 


norm 


’norm 


GF • FF 

norm norm 


(e^) = (AF„-AF } ) 

= (AF x -AF y GF EF) 
sin(A^ x /2) sm(N yWy /2) 


(16) 


N x sin(^ Y / 2) N y sin^^ / 2) 


sin (/3h cos d)EF n 


where EF norm is the normalized element factor. Note that there are 0 and (/) components 
of F nonn (determined by the component of EF norm used in Equation (16)). 


C. GRATING LOBES 

Considering a linear array along the x-axis, the array factor is a pattern that is 
symmetric about the axis of the array [9]. Therefore, the visible region is the region of the 
array factor that corresponds to - 90° <6 < +90° . The visible region is determined by the 
element spacing of the array in tenns of wavelength. One period of the array factor 
appears in the visible region when the element spacing is one-half of the wavelength. 
When the element spacing exceeds one-half wavelength, more than one period of the 
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array factor is visible and there may be more than one major lobe in the visible region. 
The different visible regions for different element spacing for an array factor with N=20 
elements is shown in Figure 6. 



Figure 6. Visible regions for different element spacing d for N = 20 element array 

(After [6]). 


These additional major lobes, which have peak intensities equal to that of the 
main lobe, are called grating lobes. The subarray spacing for a DSA is generally greater 
than one wavelength. Consider a linear DSA made up of M identical and equally-spaced 
subarrays separated by distance l x with each subarray having V-elcment equally-spaced 
half-wave dipoles with element spacing d x as shown in Figure 7. 


Subarray 1 Subarray 2 Subarray M 



2 ... N 



2 ... N 


1 2 ... N 

hi— ir nr 


hi— ir nr 

. . . 

hi— ir nr 


A 1 



dP 


h-H 

'dP 


◄- 

—► 




h 





Figure 7. Linear arrangement of equally-spaced M subarrays with each subarray 
consisting of V-clement equally-spaced half-wave dipoles. 
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The grating lobes produced by the DSA when M =5 , N = 5, l x = 1. 5 A and 
d x =1.5A are shown in Figure 8. Broadside of the array is 6 = 0°. Grating lobes are 
located at about ± 42° . 


Grating lobes 



70 80 90 


Figure 8. Grating lobes produced by linear DSA 
(M = 5, N = 5, l x = 1.5A, d x =\.5A). 

If we consider a radar array antenna with the element and subarray spacing shown 
in Figure 8, it will not be able to distinguish a target at an angle of + 42° from one at its 

broadside (0 = 0°). Therefore, the radar may be unable to associate the target with the 
correct angle, creating angular ambiguity. Another problem caused by grating lobes is the 
unwanted clutter echoes that are picked up by the grating lobes. The result is a reduction 
in the signal-to-clutter ratio which can limit target detection [10]. 

As grating lobes are generally the result of the periodicity associated with widely- 
spaced identical contiguous subarrays, one approach to suppress grating lobes is to 
implement subarrays of unequal sizes, with random locations of the subarray centers [11]. 
Another approach is to implement overlapping-subarray architecture to push the grating 
lobes away from the main beam and shape the subarray patterns in order to suppress the 
grating lobes in the subarray’s low side lobe region [12]. Another possible method to 
suppress grating lobes is to use separate transmit and receive array architectures. This 
method is based on the principle of pattern multiplication for the two-way pattern design. 
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The grating lobes in the transmit pattern can be suppressed or eliminated by placing nulls 
in the receive pattern positions coincident with the transmit grating lobes or vice versa. 

D. SUMMARY 

The fundamental theory of array antennas was presented in this chapter. The array 
pattern was obtained using the principle of pattern multiplication of the element factor, 
array factor, and ground plane factor. A two-dimensional array above a ground plane was 
considered. Grating lobes are produced when either the element spacing or subarray 
spacing is more than one-half wavelength. In most radar applications, it is undesirable to 
have grating lobes as they can cause angular ambiguities and excess background clutter. 
The subarray spacing for a DSA is generally of multiple wavelengths, and thus, grating 
lobes will exist. Therefore, it is of interest to suppress these grating lobes in DSA design. 
In the next chapter, the radar range equation for DSA, directivity and gain are introduced. 
The approach to suppress undesired grating lobes using the principle of pattern 
multiplication is also discussed. 
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III. DISTRIBUTED PHASED ARRAYS 


The principle of digital beamforming, the radar range equation for DSA, 
directivity, and gain are presented in this chapter. Undesired grating lobes produced by 
widely-spaced subarrays can be suppressed by using the principle of pattern 
multiplication. This approach is discussed in this chapter. 

A. BEAMFORMING IN PHASED ARRAYS 

Electronic scanning or beamforming in phased arrays is achieved by varying the 
amplitude and phase between antenna elements. A conventional phased array usually 
employs an analog beamforming network consisting of microwave transmission lines and 
power dividers. These are physically large and heavy if there are a large number of 
antenna elements, especially for ground or shipboard phased arrays. With the advent of 
digital processing, modern phased arrays may use digital beamforming networks. A 
phased array that employs a digital beamformer is called digital phased array. A linear 
array with digital beamforming on receive is shown in Figure 9. 


z 



Figure 9. Digital array beamforming on receive (After [13]). 
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The array factor for the linear array in Figure 9 in the x-z plane ((/> = 0° ) is given 
by 

N N 

AF{0,(j )) = ^ w n exp[ jfid n sin^cos^] = A n e^ ln exp[ jj3d n sin#cos^] (17) 
n =1 n =1 

where N is the number of elements, w n - A n e ja ' 1 is a complex weight added by the 
processor, and d n is given by Equation (1). 

If quadrature demodulators are used at each element, they provide spatial samples 
of the exponential factor in Equation (17): 

/„ = cos[j/3d n sin6*cos^] (18) 

Q n = sin[jj3d n sin6*cos^]. (19) 

Therefore, to obtain a response equivalent to a beam scanned in the direction (t? v ,^ v ), the 
weights must be given by 

w n = exp [- j/3d n sin 0 S cos (p s ]. (20) 

B. RADAR RANGE EQUATION FOR DISTRIBUTED SUBARRAYS 

The radar range equation (RRE) for distributed radar systems is derived in [4] and 
also summarized in [14]. Similarly, we can extend the RRE derivation to radars 
employing DSAs if we consider the general geometry of a DSA system as shown in 
Figure 10. 

The following assumptions are made to simplify the derivation: 

• The effects due to multipath are neglected. 

• The mutual coupling between subarrays and subarray elements is 
neglected. 

• The signals are of one polarization. 

• The case of monostatic scattering from the target applies. 
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Figure 10. General geometry of distributed subarray system (After [14]). 


If the target is in the far-field with respect to the DSA, the time-averaged (denoted 
by an over-bar) scattered power density from the target at receive subarray m when 
subarray n is transmitting is given by [ 1 ] 


W = 

rr mn 


f P t G 0 \ 

l n n 

f a ) 

V 4 ^ R n y 

v 4ftR m y 


( 21 ) 


where P t is the time-averaged power transmitted by subarray n, and G a is the antenna 

gain of subarray n. Note that this includes the element and ground plane factor for 
subarray n. The parameter (7 is the RCS of the target (assuming monostatic scattering). 

The amplitude of the peak phasor electric field at subarray m is obtained from the 
scattered power density 


W = 

rr mn 


2 % 


( 22 ) 


The phase of the electric field at receive subarray m is determined by the path length 
( R n +R m ). Therefore, the scattered electric field received by subarray m due to 
transmitting subarray n is 
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(23) 


E s = 

^mn 


2, lo P t „ 


<J 


(4n) 


2 d2 d2 

lx n A m 


exp{j +7?,„) + ^ f 


where y/ 1 is the phase added to transmitting subarray n to scan and focus the collective 
array beam. 

Hence, for N s transmitting subarrays, the total scattered electric field at the 
receiving subarray m is given by 


Ah 


Ah 


E s y E s y 


n =1 




2 *Io P tG 0 a 


(4;r) 2 RlRl 


exp {j\-P(Rn+Rm) + Vt n 


(24) 


The total received power can be obtained by multiplying the incident power 
density by the effective area of the receive subarray m, A p . Since the relationship 

between gain and effective area [9] is 



a, 


4 n 


(25) 


we can write the total received power as 


- 

p _ 

r ~ ~ . \3 


(4 k) 


N s 

Z 

n -1 


P t G 0 a 

L n u n 


d2 d2 

a /z A m 


exp {j[~0{ R n + ) + ¥t n + ¥r m ] 


(26) 


where y /is the phase added to the receive channel to focus and scan the beam. 

' m 

Since the case of monostatic scattering is considered and we can assume that 
equal phase focusing is used to achieve coherence on the target, we can make the 
following substitutions: 

• For the case of monostatic scattering: R n = R m = R a . 

• If the beam is focused on target: y/ tn +¥ r = P{ R n +R m) = 2p R 0 • 
For a single receive subarray m, the total received power can be simplified to 
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P r,„ = 


G 0 ^(7 


(4^) 3 R A 0 


AT 


'LJ p >. G o. 


n =1 


(27) 


Therefore, if N s transmit subarrays and M s receive subarrays are used, the total power 
received by the DSA is given by 

\2 

(28) 




M s 

P r = ^ P r = -i / 

Lu r m 3 4 Z_j 

m=l 14^-j A m =l 




Nr 


^K G o 

n= 1 


C. DIRECTIVITY AND GAIN 

The directivity of an array antenna is defined as the ratio of the radiation intensity 
in a certain direction to the average radiation intensity [9]. The maximum directivity is 
given as 


D(04) = 27\F mrm (e,i,f 


where Q A is the beam solid angle [9] and is defined as 


(29) 


»27i 2 

L \ F norm( 0 ’&)\ sin OdOd<j> 


(30) 


and F norm is the normalized pattern factor of the total array. Note that the nonnalized 

subarray pattern for each subarray is of the fonn of Equation (16). Substituting (30) into 
(29), we get 

D ^) = ^—, =7 --|7w,„(M 2 - (31) 

lo \o\ F norm sin dcWch/) 

From Equation (31), we observe that the directivity is determined by the total 
radiation pattern of the antenna. The power gain measures how efficient the antenna array 
is in transfonning input power to radiated power and is related to directivity by 

G(ej) = E D{ej) (32) 
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where £ is the efficiency of the antenna array. Since the antenna effective aperture A e 
[9] is related to beam solid angle by 

A 2 

(33) 

we can write the power gain of the antenna array as 

G(0,*) = ^(0,«i )| 2 . (34) 

D. PATTERN MULTIPLICATION 

When a DSA is made up of identical subarrays periodically arranged with equal 
spacing between subarrays, the total pattern of the DSA is simply the multiplication of 
the subarray pattern and the configuration pattern determined by the arrangement of the 
subarrays [5], For a planar DSA in the x-y plane composed of M x by M y uniformly 

excited subarrays with spacing l x and l y , the nonnalized pattern factor of the total array 
is 


Fnonn (e,4>) 

^norm s (9,0)xAF c (O,0) 


(35) 


where F norms ( 0, <f>) is the normalized pattern factor of the subarray from (16) and 


AF r = AF r AF v 

l x c y c 


sin(M,4/2)sin(M v ^/2) 

sin ( 4 / 2 ) sin^ 3 ? / 2 ) 


(36) 


where 

E, x = pi x (sin 6 cos (j) - sin G s cos </> s ) 

(37) 

£ y = pi y (sin 9 sin ^ - sin 6 S sin </> s ). 


Since the geometric arrangement of the subarrays is uniform, the grating lobes of 
the total array pattern can be predicted. The angular directions of the grating lobes [10] in 
the DSA configuration pattern are located at 
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(38) 


sin 6 cos (f) - sin 0 S cos <f> s = ± — p 
and 

X 

sin 6 sin (j) - sin 6 S sin <j) s = ±—p (39) 

l y 


where p,q - 0,1, 2,<x>. Note that grating lobe locations do not depend on the number 
of elements in the subarray. 

Similarly, the null locations of the uniform subarray pattern can be found at 


and 


sin 6 cos ^ 


sin 6 sin (j) 


sin 0 S cos <f> s = ± 


sin 6 s sin </> s = ± 


k 

N x d x 

A 

NyCi y 


P 


q 


(40) 


(41) 


except for p, q integer multiples of N x , N y , respectively. 

From (40) and (41), it is observed that the locations of nulls depend on the 
number of elements in the subarray. Observing (38) and (39), we see that the condition in 
which the grating lobes of the DSA configuration pattern coincide with the nulls of the 
subarray pattern is when l x / d x = N x and l y / d y = N y . However, this is the condition 

for contiguous subarray [12]. In order to suppress some ratio of the grating lobes in the 
DSA configuration pattern and still have widely-spaced subarrays congruent to DSA 
design, the subarray-to-element spacing ratio and corresponding number of elements per 
subarray need to be correctly chosen. A simple program written using MATLAB can be 
used to plot out the grating lobe and null locations in the direction cosine space for a 
selected number of elements in a subarray, the element spacing, and subarray spacing to 
help visualize the suppression of grating lobes of the configuration pattern using 
coincident nulls in the subarray pattern. The MATLAB codes for the program are given 
in Appendix B. Using an example of a transmit array with the configuration 
N x = N y =5, d x = d y = 0.51, l x = l y = 5A , the grating lobe and null locations in 
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direction cosine ( u , v ) space are shown in Figure 11. In this example, the scanning angles 
are set to 0 S =0° and <j> s = 0°. From Figure 11, we observe that the even-numbered 

grating lobes caused by the subarray spacing are suppressed by coincident nulls from the 
subarray pattern. Other grating lobes which do not have nulls that are coincident in the u, 
v space are not suppressed. 



0 

Nulls 

+ 

Grating Lobes 


Main Beam 


Figure 11. Grating lobes and null locations for transmit array with 

N x = N y =5, d x = d y = 0.5 A, l x = l y = 5 A ., 0 S = <f> s = 0°. 


By applying the principle of pattern multiplication, the normalized two-way 
pattern is defined as 

^norm 2 way ( $■> ft) ~ Fnorm Tx x ^norm ^ (42) 

where P non n Tx ( ^d/)) is the normalized total pattern of transmit array and P mrmRx ( 0,(j)) 
is the normalized total pattern of receive array. Using the two-way pattern multiplication 
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approach, we see that the remaining grating lobes can be suppressed with nulls from the 
receive array. If the receive array is configured as a contiguous array such that 
N x = N y = 10, d x = d y = 0.5/1, l x = l y = 52, the grating lobe and null locations will be 

coincident as shown in Figure 12. 



□ 

Nulls 

X 

Grating Lobes 

+ 

Main Beam 


Figure 12. Grating lobes and null locations for receive array with 

N x = Ny =10, d x = d y = 0.52, l x = l y = 52 , 0 S = <f> s = 0°. 


Overlapping the plots for the transmit and receive arrays, it is observed that the 
locations of correctly placed nulls of the receive array are coincident with the remaining 
grating lobes from the transmit array as shown in Figure 13. Therefore, the remaining 
grating lobes can be suppressed. 
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Overlap of Grating Lobes & Nulls for XMTR & RCVR Arrays. 6 S = 0°, <J> S = 0° 



0 

XMTR Nulls 

+ 

XMTR Grating Lobes 

□ 

RCVR Nulls 

X 

RCVR Grating Lobes 

+ 

Main Beam 


Figure 13. Overlap of the grating lobes and null locations for transmit and receive 

arrays. 


E. SUMMARY 

In this chapter, the principles of digital beamforming, directivity, and gain for 
phased arrays were presented, and the RRE for radar systems employing DSAs was also 
examined. It is understood that undesired grating lobes produced by widely-spaced 
subarrays can be suppressed using the principle of pattern multiplication by intentional 
placement of nulls coincident to grating lobe locations. To visualize the placement of 
grating lobes and nulls for transmit and receive array patterns, their locations are plotted 
in direction cosine space using a simple program developed in MATLAB. This forms the 
basis of two-way pattern synthesis in deciding the optimum element spacing, subarray 
spacing, and number of elements in each subarray, for the transmit and receive arrays. 
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In the next chapter, a simulation tool developed to perfonn the two-way pattern 
calculations for DSAs is described. The results of using the tool to simulate the two-way 
patterns from different DSA configurations is presented and analyzed as well. 
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IV. TWO-WAY PATTERN DESIGN FOR DSA 


The simulation tool developed to perform the two-way pattern calculation of 
transmit and receive DSAs is described in this chapter. The analysis and simulation 
results obtained using the simulation tool for different DSA configurations are also 
presented. 

A. SIMULATION TOOL FOR TWO-WAY PATTERN 

The Two-way Antenna Pattern Calculation for Distributed Subarrays program 
is a simulation tool developed and implemented in MATLAB to perform the calculation 
and display of transmit, receive, and two-way patterns for user configured DSAs. The 
graphical user interface (GUI) of the program menu is shown in Figure 14. 



Figure 14. GUI of program main menu. 

Note that the array is in the x-z plane with the y-axis normal to the aperture as 
defined by the coordinate system shown in Figure 15. The ground plane is in the x-z 
plane and the array is considered to be placed at a height above the ground plane (in the 
positive y direction). 
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Figure 15. Coordinate system with planar array in the x-z plane orientation and pattern 

angles 0 and (j> (After [8]). 

To run the program GUI, the user will need to run MATLAB and set the path or 
directory to where the program and its subroutine codes are stored. The program 
directory should contain the following MATLAB files: 

• two_way_pattern.m 

• two_way _pattern.fig 

• helpfig.m 

• helpfig.fig 

• dsaplot.m 

• caf_hclip2.m 

• caf_sdip2.m 

• caf_iso2.m 

• getamplitudes.m 

• bayliss.m 

• cosine, m 
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• tayl.m 

• compute _gain.rn 

• gausq20.m 

The MATLAB codes and description of the program files are given in Appendix 
C. After setting the path or directory, type in two_way_pattern at the MATLAB 
command line and the program GUI will be executed. 

1. DSA Parameter Inputs 

The program GUI provides a convenient way to key in the DSA parameters for 
different transmit/receive configurations. The required inputs are classified into basic 
parameters, pattern plot range settings, scan angle settings, transmit DSA configuration 
settings, and receive DSA configuration settings. 

a. Basic Parameters 

The basic parameters that the user needs to enter in the program GUI are 
frequency (in MHz), element type, dipole direction, and the height of the elements above 
the ground plane. The wavelength in meters is automatically calculated by the program 
when the array’s frequency information is entered. The user is allowed to select the 
element types: half-wave dipole, short dipole, or isotropic element from a pull-down list. 
When the half-wave dipole or short dipole is selected, the program allows the user to 
choose either v-direeled (parallel) or e-direeled (collinear) element orientations. The user 
is able to set the height of the array elements above a ground plane in terms of 
wavelengths. 


b. Pattern Plot Range and Step Size 

The program computes the DSA pattern plot over the maximum range 

0° < 6 < 180° and 0° < </> < 180°. However, the user can configure that range of angles 

for the pattern calculation under the GUI panel labeled Settings for Pattern Plot. The user 

is allowed to change the start angle, stop angle, and step size for 6 and (j> in degrees. The 

range and step size settings are applied to the transmit DSA, receive DSA, and two-way 

pattern calculation. If the start and stop values for 6 are identical, then a pattern cut for 
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that value of 0 is generated and displayed. Similarly, if the start and stop values for <j> 
are identical, then a pattern cut for that value of tj) is generated and displayed. If a range 
of both 6 and f is given, then a three-dimensional mesh plot is generated. The step size 
setting allows the user to enter an appropriate step interval for 6 and (/> angles. A small 
step size gives fine resolution data plots but requires more computation time, while a 
large step size requires less time to compute but produces coarse data plots. The default 
value for the step size is 3°. 

c. Beam Scanning Angles 

The array scan angle refers to the steering of the main beam by phase 
control. The user is able to set the scan angles for both 9 and ^ independently for the 
range 0 ° <0 < 180° and 0 ° <(j) <180° . The default scan angles are set to 90° . 

d. Transmit/Receive DSA Configuration Settings 

Two similar array GUI panels in the main menu allow the transmit and 
receive DSAs to be configured independently. Within each panel, the user is allowed to 
configure the number of subarrays (M X ,M ,) and spacing (l x ,L) in both v and z planes 
as well as select the desired amplitude distribution. The default amplitude distribution is 
uniform, but a pull-down list box allows the user to choose from taylor, cosine, bayliss, 
and triangular amplitude distributions. In the case of taylor, cosine, and bayliss 
amplitude distributions, the user is able to select the arbitrary sidelobe level (SLL). The 
number of elements (N x , N z ) and element spacing ( d x , d z ) in both v and z planes that 

form a subarray can be configured. Similarly, the amplitude distribution for the elements 
in the subarray can also be set. 

2. Two-Way DSA Pattern Calculation 

After the parameters for the transmit and receive DSAs are set, the user can select 
the Calculate Pattern button on the program GUI to initiate the computation of the array 
patterns of the transmit and receive DSAs as well as the two-way pattern. The process 


30 



flow diagram to compute and plot the array patterns is shown in Figure 16. The flow 
diagram of the sub-process to calculate the DSA one-way array pattern for half-wave 
dipoles is shown in Figure 17. 



Figure 16. Flow diagram for calculating two-way pattern of DSAs. 
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Figure 17. Flow diagram of function caf_hdip2.m to calculate DSA one-way pattern 

for half-wave dipole elements. 
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3. 


Output Plots for DSA Patterns 


When the start and stop values for 0 are identical, the program generates and 
displays the rectangular plots for the normalized transmit DSA, receive DSA, and two- 
way pattern cut for that value of 0 in both 0 and (j) components. Similarly, if the start 
and stop values for <j) are identical, the program generates and displays the rectangular 
plots for the normalized transmit DSA, receive DSA, and two-way pattern cut for that 
value of (j) in both 0 and tj) components. For the configuration settings from Figure 14, 
the rectangular plots for the normalized transmit DSA, receive DSA, and two-way pattern 
cut for (j) = 90° in 0 component are illustrated in Figures 18, 19 and 20 respectively. 
When a range of both 0 and (j) are given, the program produces normalized three- 
dimensional mesh plots for the transmit, receive, and two-way patterns of the DSAs in 
both 0 and (j) components. Using similar configuration settings from Figure 14, but 
setting the pattern plot range for both 6 and tj) from 0° to 180°and the step size 1°, we 
get the output mesh plot for the two-way pattern for the 0 component in Figure 21. The 
mesh plot is given in direction cosine (u, w ) space. 



Figure 18. Normalized transmit DSA pattern for 6 component, tj) = 90°pattern cut. 
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r igure 19. 


Normalized receive DSA pattern for 0 component, tj) = 90° pattern cut. 



Figure 20. Normalized two-way pattern for 6 component, tj) = 90° pattern cut. 
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Two-way Pattern 



Figure 21. Example of mesh plot of normalized two-way pattern for 6 component. 

4. Gain Calculation 

The user can select the Calculate Gain button on the program GUI to initiate the 
gain calculation of the transmit and receive DSAs as well as the two-way pattern. The 
computation of gain is performed based on the parameters configured on the program 
GUI. Gain calculation is performed numerically using Ng and intervals of Gaussian 

quadrature integration, with 20 points per interval, over the range 0° <0 <\ 80° and 
0° < ^ < 180° . The integration constants are loaded from MATLAB file gausq20.m. 
The user is required to set the number of integration intervals for the range of 6 and (j) in 
a user input dialog box after the Calculate Gain button is selected. The user input dialog 
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box is shown in Figure 22. The user should enter positive integer values only. The total 
number of integration points is the product of the number of integration intervals 
multiplied by 20. For example, if the integration interval for 6 is set to four, the total 
number of integration points is 80 over the range 0° <6 <\ 80°. If there are also 
four intervals in (j), the total number of integration points is 80x80 = 6400. A large 
integration interval produces more accurate gain results but increases the computation 
time. 



Figure 22. User input dialog box to set integration intervals for 6 and tj> for gain 

calculation. 

The process flow diagram for two-way gain calculation is shown in Figure 23. 
The program calls a sub-process from MATLAB function file compute_gain.m to 
calculate the one-way gain of the transmit DSA and receive DSA. The flow diagram of 
the sub-process for one-way gain calculation is shown in Figure 24. 

After the computation is completed, power gain results for the transmit DSA, 
receive DSA, and two-way pattern are displayed in both numeric and decibel format. The 
power gain results computed for the configuration settings in Figure 14 are shown in 
Figure 25. An integration interval value six is used for both 0 and (j). 
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Figure 23. Flow diagram for calculating two-way power gain of DSAs. 
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Figure 24. Flow diagram of function compute_gain.m to calculate power gain for a 

DSA. 
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Figure 25. Example of power gain results dialog window. 

5. Pattern and Configuration Data Files 

After the computation of the array patterns, the results in complex form and the 
configuration settings are auto-saved to MATLAB binary data format files 
dsapattern.mat and dsaconfig.dat, respectively, in the program directory. These fdes 
allow the user to extract the pattern data and configuration settings easily in MATLAB. 
The files are overwritten when a new pattern calculation process is initiated. 

B. DSA SIMULATIONS AND RESULTS 

To illustrate the suppression of grating lobes by placement of coincident nulls 
using the principle of pattern multiplication, the example of transmit and receive array 
configuration settings from Figure 26 is used. An isotropic element is selected as the 
antenna element type, and the transmit and receive arrays are placed at a height of 0.25 k 
above a ground plane. Uniform amplitude distribution is used for both transmit and 
receive arrays with beam scan angles 0 S = tf s = 90° (i.e., no beam scan). The transmit 

antenna consists of DSAs, and the separate receive antenna is configured to be 
contiguous designed with an array pattern that has nulls in the direction of the transmit 
array’s grating lobes. 
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Figure 26. Configuration settings for transmit and receive DSAs to show suppression 

of grating lobes. 

Using the configuration N x = A, = 5, d x = d 7 = 0.5 A, M x =M„=5,l x =L=5A 
for the transmit array, we have l x / d x = 2 N x and L / d 7 = 2N Z . Using the simulation 
tool, we obtain the normalized transmit array pattern for the 6 component with pattern 
angle range 0° < 6 <180° , (f> = 90° cut shown in Figure 27. From Figure 27, we observe 
that the even-numbered grating lobes due to subarray spacing at 6 pattern angles 36.8°, 
66.4°, 113.6° and 143.2° are suppressed by coincident nulls from the subarray pattern. 
However, the grating lobes at 0 pattern angles 11°, 53.2°, 78.7°, 101.3°, 126.8° and 169° 
are not suppressed. The peak SLL due to grating lobes is observed to be -3.7 dB. From 
Figure 27, the beamwidth between first nulls (BWFN) of the transmit array is 
approximately 4.6°. The half-power beamwidth (HPBW) of the transmit array is 
approximately 2°. The three-dimensional mesh plot and contour plot for the pattern 
angles 0° < 0 < 180° and 0° < (f) <180° of the transmit array pattern for the 0 component 
are shown in Figures 28 and 29, respectively. The step size for both 0 and (/) of the 
pattern plots are set to 0.3°. From the mesh and contour plots, we are able to observe the 

grating lobe locations in direction cosine ( u , w) space. 
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Figure 27. Normalized transmit array pattern for the 6 component, (j) = 90° cut with 
configuration settings N x = N? = 5, d x =d z =0.5A, M x =M_ = 5, l x =L = 52. 


Transmit DSA Pattern 



U = sin0*cos<(> 


Figure 28. Normalized mesh plot of transmit array pattern for the 0 component with 
configuration settings N x =N Z =5, d x =r/_ =0.5A, M x =M_ = 5, l x =L= 52. 


42 

































































































Transmit DSA Pattern 
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Figure 29. Normalized contour plot of transmit array pattern for the 9 component with 
configuration settings N x = N z =5, d x =cL =0.5/1, M x =M_ =5, l x -L =5/1. 

To place nulls at the remaining grating lobe locations of the transmit array, we 
can apply a contiguous subarray configuration for the receive array with the settings 
given in Figure 26. The normalized receive array pattern for the 9 component with 
pattern angle range 0° <9 < 180°, (/) = 90° cut is shown in Figure 30. Since the receive 
array is contiguous with element spacing 0.5/1, no grating lobes exist in the visible 
region. The peak SLL is observed to be -13.3 dB. The BWFN of the receive array is 
approximately 4.6°, and FIPBW of the receive array is approximately 2°. The three- 
dimensional mesh plot and contour plot for the pattern angles 0° < 9 < 180° and 
0 ° <(j) <180° of the receive array pattern for the 9 component are shown in Figures 31 
and 32, respectively. The step size for both 9 and (/) of the pattern plots are set to 0.3°. 
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Figure 30. Normalized receive array pattern for the 6 component, (f) = 90° cut with 
configuration settings N x =/V_ =10 ,d x =d z =0.5/1, M x -M„ =5 ,l x —l 7 =5/1. 


Receive DSA Pattern 



U = sin0*cos4> 

Figure 31. Normalized mesh plot of receive array pattern for the 6 component with 
configuration settings N x = /V_ =10, d x = d~ =0.5/1, M x -M„ =5, l x -L =5/1. 
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Figure 32. Normalized contour plot of transmit array pattern for the 9 component with 

configuration settings N x =N 7 =10, d x =cL =0.5/1, M x =M„ = 5, l x = l z - 5A. 

The normalized two-way pattern of the 6 component, (j) = 90° pattern cut, for the 
configuration settings given in Figure 26 is shown in Figure 33. ft is observed that the 
grating lobes at 0 pattern angles 11°, 53.2°, 78.7°, 101.3°, 126.8° and 169° of the 
transmit array have been suppressed. The peak SLL of the two-way pattern is -25.5 dB 
due to the subarray sidelobes at 6 pattern angles 86.7° and 93.3°. The grating lobes at 
0 = 78.7°, 101.3° are suppressed to -29 dB. The improvement factor in grating lobe 
suppression is 25.3 dB. The BWFN of the two-way pattern is approximately 4.6°, and 
HPBW of the two-way pattern is approximately 1.5°. The three-dimensional mesh plot 
and contour plot for the pattern angles 0 ° <9 < 180° and 0° < (/) <180° of the two-way 
pattern for the 9 component are shown in Figures 34 and 35, respectively. The step size 
for both 9 and (/) of the pattern plots are set to 0.3°. 
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Figure 33. Normalized two-way pattern for the 6 component, (f> = 90° cut of a 
transmit DSA with contiguous receive subarrays configuration. 


Two-way Pattern 
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Figure 34. Normalized mesh plot of two-way pattern for the 6 component in direction 

cosine space. 
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Figure 35. Normalized contour plot of two-way pattern for the 6 component in 

direction cosine space. 


Two-way Pattern 



Using an integration interval of ten for both 0 and (/>, we get the power gain 
results of the transmit array, receive array, as well as the two-way power gain shown in 
Figure 36. The simulation results of the DSA configuration settings in Figure 26 are 
summarized in Table 1. 


S3 Gain Results 


T ransmit Array: 

Numeric gain = 1955.4256, Gain in dB = 32.9124 dB 
Receive Array: 

„ Numeric gain = 7797.4208, Gain in dB = 38.9195 dB 

T wo-way numeric gain = 15247276.1548 
Two-way Gain in dB = 71.8319 dB 


OK 


Figure 36. Power gain results of transmit DSA with 
contiguous receive subarrays configuration. 
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Table 1. Summary of simulation results of transmit DSA with contiguous receive 

subarrays configuration. 



Peak SLL 

BWFN 

HPBW 

Power Gain 

Transmit Array 

-3.7 dB 

4.6° 

2° 

32.9 dB 

Receive Array 

-13.3 dB 

4.6° 

2° 

38.9 dB 

Two-way Pattern 

-25.5 dB 

4.6° 

1.5° 

71.8 dB 


From the simulation results, it is observed that two-way pattern multiplication 
approach is effective in suppression of grating lobes. In this example, a grating lobe 
suppression of at least 25 dB can be achieved. 

In the next example, both transmit and receive antennas are configured as widely 
spaced, thinned arrays [15] to achieve a narrow scanned beam for very good angular 
resolution without the cost and space required of a fully filled array. Chebyshev 
illumination is applied to achieve low sidelobe patterns. The transmit array and receive 
array configuration settings is given in Figure 37. 



Figure 37. Configuration settings for thinned transmit and receive arrays. 
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The transmit array is made up of 16x16 subarrays with 3 A spacing. The receive 
array is made of 16x16 subarrays with 4A spacing. Each subarray consists of 4x4 z- 
directed half-wave dipole elements placed at a height of 0.25 A above a ground plane. 
Taylor amplitude distribution with 40 dB SLL and n - 5 is applied to both transmit and 
receive array configurations. Beam scan angles 6 S = (j s — 90° are considered. The 
normalized transmit array pattern of the 16x16 subarrays for the 6 component with 
pattern angle range 0° < 0 < 180° , (jr= 90° cut is shown in Figure 38. The peak SLL due 
to grating lobes at 6 = 70.6°, 109.4° is observed to be -7.9 dB. The BWFN and HPBW 
of the transmit array pattern is approximately 4.4° and 1.5°, respectively. The three- 
dimensional mesh plot and contour plot for the pattern angles 0° < 0 < \ 80° and 
0 ° <(j) <180° of the transmit array pattern for the 0 component are shown in Figures 39 
and 40, respectively. A step size of 0.1° is used for both 6 and (j) of the pattern plots. 



Figure 38. Normalized transmit array pattern of 16x16 subarrays with spacing 3 A for 

the 6 component, ^ = 90° cut. 
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Transmit DSA Pattern 



U = sir0*cos(|> 


Figure 39. Normalized mesh plot of transmit array pattern of 16x16 subarrays with 

spacing 3 A for the 0 component. 


Transmit DSA Pattern 



U = sirfl*cos4> 


Figure 40. Normalized contour plot of transmit array pattern of 16x16 subarrays with 

spacing 3 A for the 6 component. 
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The normalized receive array pattern of the 16x16 subarrays for the 6 
component with pattern angle range 0 ° <6 < 180° , (f) = 90° cut is shown in Figure 41. 
The peak SLL due to grating lobes at 6 = 75.6°, 104.4° is observed to be -4.1 dB. Notice 
that the grating lobe locations of the transmit and receive arrays are not coincident due to 
their different subarray spacing. The BWFN and HPBW of the receive array pattern is 
approximately 3.2° and 1.1°, respectively. The three-dimensional mesh plot and contour 
plot for the pattern angles 0° < 6 < 180° and 0° <(/)< 180° of the receive array pattern for 
the 6 component are shown in Figures 42 and 43, respectively. A step size of 0.1° is 
used for both 6 and (/) of the pattern plots. 


Receive DSA Pattern. 0=90° cut 



Figure 41. Normalized receive array pattern of 16x16 subarrays with spacing 4/L for 

the 6 component, (j) = 9 0° cut. 
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Receive DSA Pattern 



W = cos8 


U = sin8*cos4> 


Figure 42. Normalized mesh plot of receive array pattern of 16x16 subarrays with 4A 

for the 6 component. 
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Figure 43. Normalized contour plot of receive array pattern of 16x16 subarrays with 

4/1 for the 0 component. 
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The normalized two-way pattern of the 6 component, (f) = 90° pattern cut, is 
shown in Figure 44. It is observed that the grating lobes for both transmit and receive 
arrays have been suppressed since their locations are not coincident, and both transmit 
and receive array patterns have low SLLs (exclusive of the grating lobes). The peak SLL 
of the two-way pattern is -49.7 dB at 0 = 76.1°, 103.9°. The BWFN of the two-way 
pattern is approximately 4.4°, and HPBW of the two-way pattern is approximately 0.9°. 
The three-dimensional mesh plot and contour plot for the pattern angles 0° < 0 < 180° 
and 0 ° <(j) <180° of the two-way pattern for the 0 component are shown in Figures 45 
and 46, respectively. The step size for both 6 and (f> of the pattern plots are set to 0.1°. 


Two-way Pattern, ^=90° cut 



Figure 44. Normalized two-way pattern of thinned transmit and receive arrays for the 

0 component, ^ = 90° cut. 
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Two-way Pattern 
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Figure 45. Normalized mesh plot of two-way pattern of thinned transmit and receive 

arrays for the 6 component. 
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Figure 46. Normalized contour plot of two-way pattern of thinned transmit and receive 

arrays for the 6 component. 


Two-way Pattern 
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Using an integration interval of 18 for both 6 and (j), we get the power gain 
results of the transmit array, receive array, as well as the two-way power gain shown in 
Figure 47. The simulation results are summarized in Table 2. 



Figure 47. Power gain results using thinned transmit and receive array configurations. 


Table 2. Summary of simulation results for thinned transmit and receive array 

configurations. 



Peak SLL 

BWFN 

HPBW 

Power Gain 

Transmit Array 

-7.9 dB 

4.4° 

1.5° 

39 dB 

Receive Array 

-4.1 dB 

3.2° 

1.1° 

39 dB 

Two-way Pattern 

-49.7 dB 

4.4° 

o 

Id 

o 

78 dB 


The transmit array pattern, receive array pattern and two-way pattern when the 
main beam is scanned to 0 S =100° are shown in Figures 48, 49 and 50, respectively. It is 
observed that the peak SLL, BWFN and HPBW are typically unchanged for small beam 
scan angles. 
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Transmit DSA Pattern. 4>=90° cut 



Figure 48. Normalized transmit array pattern of 16x16 subarrays with spacing 3 A 

when scanned to 0 S = 100° . 


Receive DSA Pattern, <J)=90° cut 



Figure 49. Normalized received array pattern of 16x16 subarrays with spacing A A 

when scanned to 0 S = 100°. 
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Two-way Pattern, 4>=90° cut 



Figure 50. Normalized two-way pattern when scanned to O s = 100° . 

The transmit array pattern, receive array pattern and two-way pattern when the 
main beam is scanned to 0 S = 140° are shown in Figures 51, 52 and 53, respectively. For 

large scan angles, beam broadening effect is more evident, and the peak SLL is increased. 
The peak SLL for the two-way pattern is -40.2 dB when the beam is scanned to 
0 S = 140°. 


Transmit DSA Pattern. 4>=90° cut 



Figure 51. Normalized transmit array pattern of 16 x 16 subarrays with spacing 3 A 

when scanned to 0 S = 140° . 
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Receive DSA Pattern, $=90° cut 



0 20 40 60 80 100 120 140 160 180 

Pattern Angle, 6 (deg) 


Figure 52. Normalized received array pattern of 16x16 subarrays with spacing 4A 

when scanned to O s = 140° . 


Two-way Pattern. <t>=90° cut 



Figure 53. Normalized two-way pattern when scanned to O s = 140°. 

The two-way patterns of beams scanned at increments of A<9 ? =5° for 
75° < O s < 105° are illustrated in Figure 54. 
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Two-way Pattern, 4>=90° cut 



Figure 54. Normalized two-way pattern of multiple beam scanning for 75° < Q s < 105° 

with A0 S =5° 

We see from the simulation results that using separate transmit and receive array 
patterns with non-coincident grating lobe locations and low sidelobe illumination, the 
grating lobes can be effectively suppressed. Although beam broadening usually 
accompanies low sidelobe illuminations, a narrow beamwidth of the two-way pattern can 
still be achieved with wide subarray spacing. 

C. SUMMARY 

A simulation tool was developed and implemented in MATLAB to perform the 
two-way pattern calculation of transmit and receive DSAs. The features of the simulation 
tool were described in this chapter. The simulation tool was also able to perform gain 
calculation for the transmit and receive DSAs, as well as compute the two-way gain. 

Using the simulation tool, we demonstrated the effectiveness of the two-way 
pattern multiplication approach to suppress undesired grating lobes by placement of 
subarray nulls at the grating lobe locations. In this chapter, we also showed that low SLL 
and narrow FIPBW of the two-way antenna pattern can be achieved with thinned transmit 
and receive arrays consisting of widely-spaced subarrays with non-coincident grating 
lobe locations. 
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V. SUMMARY AND RECOMMENDATIONS 


A. SUMMARY 

The primary objective of this research was to develop a simulation tool to 
investigate the behavior and effectiveness in suppressing undesired grating lobes using 
the approach of two-way antenna pattern in DSA design. The fundamental array theory 
and principle of pattern multiplication, which fonned the basis of the simulation tool 
design, were discussed. Using the principle of pattern multiplication, we can suppress 
grating lobes by placement of subarray nulls at grating lobe locations. A simple program 
was developed in MATLAB to allow the user to visualize the placement of grating lobes 
and nulls in direction cosine space for a DSA configuration in the visible region. A 
simulation tool with GUI was developed and implemented in MATLAB to perform the 
two-way antenna pattern and power gain calculations for user configured DSAs. The 
program is capable of performing two-way pattern and power gain calculations for linear 
or planar DSAs consisting of isotropic elements, half-wave dipoles or short dipoles above 
a ground plane. The program is able to present the simulation results in the 6 pattern cut, 
(j) pattern cut, or as a three-dimensional mesh plot in direction cosine space. The program 
GUI provides a convenient way for the user to tweak the design configurations very 
quickly by changing the DSA parameters. The program can serve as a useful tool for both 
students and electromagnetic professionals to detennine and study the two-way pattern 
and power gain of different transmit/receive DSA designs. 

The effectiveness of the two-way pattern multiplication approach to suppress 

undesired grating lobes by placement of subarray nulls at the grating lobe locations was 

demonstrated using the simulation tool. For a transmit DSA with subarray spacing of 5/1 

and uniform amplitude illumination, a SLL of -25.5 dB was achieved. Using the 

simulation tool, we demonstrated that low SLL and narrow HPBW of the two-way 

antenna pattern can be achieved using thinned transmit and receive arrays consisting of 

widely-spaced subarrays with non-coincident grating lobe locations. The simulation 

results showed that a SLL of -49.7 dB and HPBW of 0.9° were obtained for broadside 

illumination. Simulation was carried out to examine the effects on the two-way pattern 
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when different scan angles were applied. The simulation results showed that beam 
broadening was evident and peak SLL was increased for the two-way antenna pattern 
when large scan angles were applied to the main beams of the transmit and receive 
arrays. It was assumed that the main beams of the transmit and receive arrays were 
scanned to the same angle. 

B. RECOMMENDATIONS 

At this time, the simulation tool allows us to study the two-way pattern for 
equally-spaced linear or rectangular DSAs. Other common array configurations include 
triangular arrays and circular arrays. A triangular grid is known to be more efficient than 
a rectangular grid for the suppression of grating lobes because fewer elements are 
required for a given aperture size [10]. Circular arrays, in which the elements are placed 
in a circular ring, have the advantage of symmetry in azimuth, making them suitable for 
applications that require full 360° coverage, such as direction-finding, air and space 
navigation, and underground propagation [16]. It would be useful to upgrade the 
simulation tool to be able to compute the two-way pattern of triangular, circular and other 
array configurations to allow the study of their two-way patterns. 

Phase shifters are used to control the location and shape of the antenna beam. In 
practice, most phase shifters are digitally controlled devices, which allow only discrete 
values for the phase shift. Therefore, truncation or round-off error must be introduced 
onto the phase shifts, and this yields a periodic quantization error that gives rise to 
quantization lobes [10]. The current version of the simulation tool assumes continuous 
phase shift in the calculation and does not account for the effect of phase quantization. It 
is essential to include some round-off error or truncation methods for the phase shift 
computation as a future improvement to the program. 
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APPENDIX A. DIPOLES OF ARBITRARY ORIENTATION 1 


The dipole is centered at the origin as shown in Figure Al. The dipole direction is 
given by the vector 


/ = x sin 6 a cos (j) a + y sin 6 a sin (j) a + z cos 6 a . 


w„ 


The dipole direction cosines are 


^-direction cosine: u a = cos a x = sin 0 a cos (j) a 
y-direction cosine: v a = cos a y = sin 0 a sin <j) a 
z-direction cosine: w a = cos ay = cos 6 a . 


(43) 


(44) 


The length of the dipole is L, and for a half-wave dipole the current distribution is given 
by the function of the path length variable t 

lU') = I m cos(j3f), -L/2<f<L/2 (45) 


where I m is the maximum value at the feed point. For an ideal (Hertzian) dipole with 
constant current 


lU') = I m , -L/2 < £' < L/2, L«A. (46) 

The direction cosines for the far-field observation point are 

x-direction cosine: u = sin 6 cos </> 

y-direction cosine: v = sin 0 sin (j) (47) 

z-direction cosine: w = cos 6. 


1 The material for Appendix A is from “Arbitrary dipoles” unpublished notes of David C. Jenn. 
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Figure A1. Dipole of arbitrary orientation centered at the origin 
The far-Held radiation integral can be written as 




-iPVoi c - jpr 
4nr 



(48) 


where it is understood that the r component is discarded. For the ideal dipole the integral 
is 


LI 2 

I 

-L/2 


J lU')e jpe ' hl) di' = I m Lsmc 


f n j j A 

pLvl 


I m L (L«A), 


and the resulting spherical field components are 


E e =E{0, </i)»0= l^JohnL e -jP r i*e 
4 nr 

E4 =E(0, </>)•</>= ~^ T] ° ImL e ~ j P r I • j). 

• A.mr 


(49) 


(50) 
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From the transform tables: 


6 = x cos 6 cos 6 + v cos 6 sin <b - z sin 6 

„ „ ( 51 ) 

<j) = -xsm<j)+ ycostj). 

For a Hertzian dipole along the z-axis ($ a = (j) a = 0°, 7 = z): 


I • 6 = z*0 = -sin 6 
/ • ^ = £•^ = 0, 


(52) 


and the fields reduce to the familiar results 


E e = jP% I m L e -jpr sin £ 

4 nr 

Ef= 0. 

For a half-wave dipole (L = 4, / 2), the integral from (48) is 


(53) 


L/2 


| I m cos(pP)e j ^ hI ' > df = 


-L/2 


I c 


COS 


[pl'{r'l)\ +p 2 


L/2 


(54) 


J-L/2 . 


After applying some trigonometric identities, we get the final result for the far electric 
field 


E(0,t) = - 


Mje-V' C 0 S U f - 7 


n 


2nr i -(?./) 2 

For a half-wave dipole along the z-axis {o a = tj) a = 0°, / = z): 


r •/ = cos 6 = co, 

0»i = -sm0, 


(55) 


(56) 


En = 


iVo^m e 


~jpr 


2 nr 


cos —cos 9 
\ 2 

sin# 


Ej= 0. 


(57) 
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For a half-wave dipole along the y-axis (o a = <j> a = 90°, / = y ): 


Ee = 



r»I 

6»i 

(j)»i 

= sin 0 sin (j) = v , 

= cos 0 sin (j), 

= COS (j) , 

m 0 hi e ~ if)r 

cos|^sin#sin^j 

2 nr 

1-sin 2 #sin 2 <f> 

Wm e ~ Jpr 

coslysin^sin^j 

2 nr 

1-sin 2 #sin 2 (/> 


cos 6 sin (f >, 


cos (j) . 


For a half-wave dipole along the v-axis = 90 °,<j> a = 0°,/ = x) ■ 


r • I = sin 0 cos (j) - u , 
§• I = cos 6 cos $ , 
^•/ = -sin^, 


Eo = 

E,/, = 


Plo I m e 


~iP r 


2 nr 


cos I y sin 0 cos (f) 
1-sin 2 <9 cos 2 (f) 


.19(J m e 


~jpr 


2nr 


cos (y sin 6 cos </> 
1-sin 2 9 cos 2 <f> 


cos 6 cos <j>. 


sin^. 


(58) 


(59) 


(60) 


(61) 
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APPENDIX B. MATLAB CODE FOR PLOTTING GRATING 
LOBES AND NULL LOCATIONS OF ARRAY PATTERN 


% Plot_rect_nulls_and_grating_xyplane.m 

% This code computes and plots the nulls and grating lobes locations 
% for rectangular lattice of DSAs in x-y plane in the direction 
% cosine space. Comparison is made between XMTR and RCVR DSAs to 
% observe the placement of nulls at grating lobe locations to 
% suppress grating lobes. 

% Modified by Cher Hock Hin 
% Date: 8 Aug 2012 

clear 

rad=pi/180; 

%XMTR array settings 

dxt=0.5; dyt=0.5; %element spacing in x and y plane (in wavelengths) 
Nxt=5; Nyt=5; %number of elements in x and y plane 
lxt=5; lyt=5; %subarray spacing in x and y plane (in wavelengths) 
%RCVR array settings 

dxr=0.5; dyr=0.5; %element spacing in x and y plane (in wavelengths) 

Nxr=10; Nyr=10; %number of elements in x and y plane 

lxr=5; lyr=5; %subarray spacing in x and y plane (in wavelengths) 

%scan angle 

thetas=0; thetar=thetas*rad; 
phis=0; phir=phis*rad; 

us = sin (thetar)*cos(phir); vs = sin(thetar)*sin(phir); 

fmhz=100; 

f=fmhz*le6; 

wave=3e8/f; 

wdx=l/dxt; %lambda/spacing 
wdz=l/dyt; %lambda/spacing 

N=20; %sets the number of grating lobes to plot 
pt=0; lt=0; pr=0; lr=0; 
for n=-N:N; 

for m=-N:N 

Dxt=n/dxt/Nxt+us; Dyt=m/dyt/Nyt+vs; 
xlobet=n/lxt+us; ylobet=m/lyt+vs; 
lobet=sqrt(xlobet A 2+ylobet A 2); 
nullt=sqrt(Dxt A 2+Dyt A 2); 
if nullt<=l %find XMTR array nulls 
pt=pt+l; 

Xt(pt)=Dxt; Yt(pt)=Dyt; 

end 

if lobet<=l %find XMTR array grating lobes 
lt=lt+l; 

Xlobet (It)=xlobet; Ylobet(It)=ylobet; 

end 

Dxr=n/dxr/Nxr+us; Dyr=m/dyr/Nyr+vs; 
xlober=n/lxr+us; ylober=m/lyr+vs; 
lober=sqrt(xlober A 2+ylober A 2); 
nullr=sqrt(Dxr A 2+Dyr A 2); 
if nullr<=l %find RCVR array nulls 
pr=pr+l; 

Xr(pr)=Dxr; Yr(pr)=Dyr; 

end 

if lober<=l %find RCVR array grating lobes 
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lr=lr+l; 

Xlober(lr)=xlober; Ylober(lr)=ylober; 

end 

%main beam 

if m==0 & n==0, X0t=Dxt; Y0t=Dyt; end 

end 

end 

%unit circle 

t=0:360; cx=cos(t*rad); cy=sin(t*rad); 
figure (1), elf %plot for XMTR array 

plot (Xt,Yt, 'bd' ,Xlobet,Ylobet, 'm+' ,XOt,YOt, ' r* ' , cx, cy, '-r' ) 

legend( 'Nulls'Grating Lobes', 'Main Beam ','Location', 'EastOutside' ) 

axis ( [-1,1,-1,1]) 

title ({[' XMTR Array, \theta_s= ', num2str(thetas ),' A o, \phi_s = 

',num2str(phis),' A o'];... 

['N_x = ',num2str(Nxt),', d_x = ',num2str(dxt),'\lambda, ','l_x = 

',num2str(lxt),'\lambda, ',... 

'N_y = ' ,num2str(Nyt),' , d_y = ' ,num2str(dyt),' \lambda, ' , 'l_y = 

',num2str(lyt),'\lambda'];}) 
axis square 

xlabel( 'u=sin\theta cos\phi') 
ylabel ( 'v=sin\theta sin\phi') 

figure(2), elf %plot for RCVR array 

plot(Xr,Yr, 'ks' ,Xlober,Ylober, 'bx', XOt,YOt,' r*' ,cx,cy, '-r' ) 

legend ( 'Nulls ', 'Grating Lobes', 'Main Beam ', 'Location', 'EastOutside' ) 

axis([-1,1,-1,1]) 

title ({[' RCVR Array, \theta_s= ',num2str(thetas ),' A o, \phi_s= 

',num2str(phis),' A o'];... 

['N_x = ',num2str(Nxr),', d_x = ',num2str(dxr),'\lambda, ','l_x = 

',num2str(lxr), '\lambda, ' , . . . 

'N_y = ' ,num2str(Nyr),' , d_y = ' ,num2str(dyr),' \lambda, ' , 'l_y = 

',num2str(lyr),'\lambda'];}) 
axis square 

xlabel (' u=sin\theta cos\phi') 
ylabel (' v=sin\theta sin\phi') 

figure(3), elf %plot overlap of XMTR and RCVR array 

plot (Xt,Yt, 'bd' ,Xlobet,Ylobet, 'm+' ,Xr,Yr, 'ks', Xlober,Ylober, ' bx' , . .. 

XOt,YOt, 'r*' ,cx,cy, '-r') 

legend ('XMTR Nulls', 'XMTR Grating Lobes', 'RCVR Nulls', 'RCVR Grating 
Lobes' , 'Main Beam ','Location', 'EastOutside' ) 
axis([-1,1,-1,1]) 

title (['Overlap of Grating Lobes & Nulls for XMTR & RCVR Arrays, \theta_s 
',num2str(thetas) ,' A o, \phi_s= ' ,num2str(phis), ' A o ' ] ) 
axis square 

xlabel (' u=sin\theta cos\phi') 
ylabel (' v=sin\theta sin\phi') 
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APPENDIX C. MATLAB CODE FOR SIMULATION TOOL 


function dsaplot(action) 

% dsaplot.m 
% Version: 2.1 
% Author: Cher Hock Hin 
% Advisor: Professor David C. Jenn 
% Date: 13 August 2012 

% Uses the function "caf_hdip2.m", "caf_sdip2.m" and "caf_iso2.m 

% to compute XMTR and RCVR array patterns for half wave dipoles, 
% short dipoles and isotropic elements which return outputs in 
% complex form (non-dB). 

% Two-way pattern is then computed using pattern multiplication. 
% Two-way gain is also computed using method of numerical 
% integration. 

% Pattern and gain calculation is based on user configuration 
% of XMTR and RCVR planar arrays consisting of dipoles or 
% isotropic elements above a ground plane. The arrays are 
% oriented in the x-z plane; ground plane is the x-z plane; 

% y axis is normal to ground plane. 

switch (action) 

%Case for frequency 
case 'frequency' 

h_frequency = findobj(gcf,' Tag ',' frequency ') ; 

freq_str = get(h_frequency, 'String ') ; 

freq = getFreq(freq_str); 

set(h_frequency, 'String' ,num2str(freq)); 

wavelength = 3e8/freq/le6; 

h_wavelength = findobj(gcf, 'Tag','wavelength'); 
set(h_wavelength, 'String' ,num2str(wavelength)); 

%Case for element type selection 
case 'eltype' 

h_eltype = get(findobj(gcf ,' Tag ',' eltype '), 'Value ') ; 
if h_eltype == 3 

set (findobj(gcf, ' Tag' , ' dipdir '), 'Enable', 'off'); 

else 

set(findobj(gcf,' Tag' ,' dipdir' ), 'Enable' , 'on' ); 

end 

%Case for phi start angle 
case 'pstart' 

h_pstart = findobj(gcf ,' Tag ',' pstart ') ; 
pstart_str = get(h_pstart ,' String ') ; 
pstart = getPStart(pstart_str); 
set(h_pstart, 'String' ,num2str(pstart)); 

%Case for phi stop angle 
case 'pstop' 

h_pstop = findobj(gcf,' Tag ', 'pstop ') ; 
pstop_str = get(h_pstop, 'String ') ; 
pstop = getPStop(pstop_str); 
set(h_pstop,' String' ,num2str(pstop)); 

%Case for theta start angle 
case 'tstart' 
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h_tstart = findobj(gcf ,' Tag' , 'tstart ') ; 
tstart_str = get(h_tstart, 'String ') ; 
tstart = getTStart(tstart_str); 
set(h_tstart, 'String' ,num2str(tstart)); 

%Case for theta stop angle 
case 'tstop' 

h_tstop = findobj(gcf,' Tag ',' tstop ') ; 
tstop_str = get(h_tstop, 'String 1 ); 
tstop = getTStop(tstop_str); 
set(h_tstop, 'String' ,num2str(tstop)); 

%Case for thetas_slider 
case 'thetas_slider' 

h_thetas_slider = findobj(gcf,' Tagthetas_slider ') ; 
thetas = ceil(get(h_thetas_slider, 'Value ')) ; 
h_thetas = findobj(gcf,' Tag ',' thetas ') ; 
set(h_thetas, ' String 1 ,num2str(thetas)); 

case 'thetas' 

h_thetas = findobj(gcf,' Tag ',' thetas ') ; 
thetas = str2num(get(h_thetas,' String' )); 

thetas = getScanAngle(thetas); %validate range between 0 to 180 
h_thetas_slider = findobj(gcf,' Tag ',' thetas_slider ') ; 
set(h_thetas_slider,' Value ceil(thetas)); 
set(h_thetas,' String' ,num2str(thetas)); 

%Case for phis_slider 
case 'phis_slider' 

h_phis_slider = findobj(gcf ,' Tag ', 'phis_slider ') ; 
phis = ceil(get(h_phis_slider, 'Value ')) ; 
h_phis = findobj(gcf ,' Tag ',' phis ') ; 
set(h_phis, 'String' ,num2str(phis)); 

case 'phis' 

h_phis = findobj(gcf ,' Tag ',' phis '); 
phis = str2num(get(h_phis,' String ')) ; 

phis = getScanAngle(phis); %validate range between 0 to 180 
h_phis_slider = findobj(gcf ,' Tag ', 'phis_slider ') ; 
set(h_phis_slider, 'Value' ,ceil(phis)); 
set(h_phis, 'String' ,num2str(phis)); 

%Case for theta step size 
case 'delt' 

h_delt = get(findobj (gcf, 'Tag', 'delt'), 'Value'); 

if h_delt >= 3, warndlg ('Setting step size to < 1 degree will greatly 
increase computation time! ', 'Warning! ') ; end 

%Case for phi step size 
case 'delp' 

h_delp = get(findobj(gcf, 'Tag’ , 'delp' ), 'Value '); 

if h_delp >= 3, warndlg (' Setting step size to < 1 degree will greatly 
increase computation time! ', 'Warning! ') ; end 

%Case for Transmit Array parameters 
case 'Nxt' 

h_Nxt = findobj(gcf,' Tag ', 'Nxt ') ; 
nxt_str = get(h_Nxt, 'String' ); 

Nxt = getNEL(nxt_str); 
set(h_Nxt,' String' ,num2str(Nxt)); 
if Nxt == 1 
dxt = 0; 
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set(findobj(gcf, 'Tag', 'dxt' ) , 'String', num2str(dxt)); 
warndlgt'For a single element, array spacing is set to zero.',... 

' Array Spacing ',' help '); 
elseif Nxt > 1 

dxt = str2num(get(findobj(gcf, 'Tag' , 'dxt' ), 'String ')) ; 
if dxt == 0 

set(findobj(gcf, 'Tag' , 'dxt' ), 'String' ,num2str(1)); 
warndlg(’For the number of elements, array spacing cannot be 

zero .',... 

' Array Spacing ',' help '); 

end 

end 

case 'Nzt' 

h_Nzt = findobj(gcf,' Tag' , 'Nzt' ); 
nzt_str = get(h_Nzt,' String ') ; 

Nzt = getNEL(nzt_str); 
set(h_Nzt,' String' ,num2str(Nzt)); 
if Nzt == 1 
dzt = 0; 

set(findobj(gcf, 'Tag', 'dzt'), 'String', num2str (dzt)); 
warndlgt'For a single element, array spacing is set zero.',... 

' Array Spacing ',' help '); 
elseif Nzt > 1 

dzt = str2num(get(findobj(gcf ,' Tag ',' dzt '),' String ')) ; 
if dzt == 0 

set(findobj(gcf, 'Tag','dzt'),'String', num2str(1)); 
warndlgt'For the number of elements, array spacing cannot be 

zero .',... 

' Array Spacing' , 'help' ); 

end 

end 

case 'dxt' 

h_Dxt = findobj(gcf,' Tag ', 'dxt ') ; 
dxt_str = get(h_Dxt,' String' ); 
dxt = getSpacing(dxt_str); 
set(h_Dxt, 'String' ,num2str(dxt)); 
if dxt == 0 
Nxt = 1; 

set(findobj(gcf, 'Tag', 'Nxt' ), 'String' ,num2str(Nxt)); 
warndlg (' Number of array elements set to 1. ', ' Number of Array 
Elements' , 'help' ); 
end 

Nxt = str2num(get(findobj(gcf,' Tag' , 'Nxt' ),' String' )); 
if dxt > 0 

if Nxt == 1 

errordlg( 'For the spacing indicated, number of elements must be at 
least 2.' ,... 

'Number of Elements' , 'error' ); 
dxt = 0; 

set(h_Dxt, 'String' ,num2str(dxt)); 

end 

end 

case 'dzt' 

h_Dzt = findobj(gcf,' Tag ',' dzt ') ; 
dzt_str = get(h_Dzt,' String ') ; 
dzt = getSpacing(dzt_str); 
set(h_Dzt, 'String' ,num2str(dzt)); 
if dzt == 0 
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Nzt = 1; 

set(findobj(gcf, 'Tag','Nzt'),'String', num2str(Nzt)); 
warndlg (' Number of array elements set to 1.',... 

' Number of Array Elements ',' warn ') ; 


end 


Nzt = str2num(get(findobj(gcf, 'Tag','Nzt'),'String')); 
if dzt > 0 

if Nzt == 1 

errordlg( 'For the spacing indicated, number of elements must be at 
least 2 . ' , . . . 

'Number of Elements' , 'error' ); 
dzt = 0; 

set(h_Dzt, 'String' ,num2str(dzt)); 

end 

end 


case 'Mxt' 

h_Mxt = findobj(gcf,' Tag ', 'Mxt ') ; 
mxt_str = get(h_Mxt, 'String' ); 

Mxt = getNEL(mxt_str); 
set(h_Mxt,' String' ,num2str(Mxt)); 
if Mxt == 1 
Ixt = 0; 

set(findobj(gcf, 'Tag', 'lxt' ), 'String', num2str(Ixt)); 
warndlg ('For a single element, array spacing is set to zero.',... 
' Array Spacing ',' help '); 
elseif Mxt > 1 

lxt = str2num(get(findobj(gcf, 'Tag ',' lxt '),' String ')); 
if lxt == 0 

set(findobj(gcf, 'Tag' , 'lxt' ), 'String' ,num2str(1)); 
warndlg('For the number of elements, array spacing cannot be 

zero .',... 

' Array Spacing' , 'help' ); 

end 

end 


case 'Mzt' 

h_Mzt = findobj(gcf,' Tag ', 'Mzt ') ; 
mzt_str = get(h_Mzt,' String' ); 

Mzt = getNEL(mzt_str); 
set(h_Mzt, 'String' ,num2str(Mzt)); 
if Mzt == 1 
Izt = 0; 

set(findobj(gcf, 'Tag', 'Izt' ), 'String', num2str (Izt)); 

warndlg ('For a single element, array spacing is set zero.',... 

' Array Spacing ',' help ') ; 
elseif Mzt > 1 

Izt = str2num(get(findobj(gcf ,' Tag ',' Izt '),' String ')); 
if Izt == 0 

set(findobj(gcf, 'Tag', 'Izt'), 'String', num2str(1)); 

warndlg ('For the number of elements, array spacing cannot be 

zero .',... 

' Array Spacing ',' help ') ; 

end 

end 


case 'lxt' 

h_Lxt = findobj(gcf ,' Tag ',' lxt ') ; 
lxt_str = get(h_Lxt, 'String' ); 
lxt = getSpacing(lxt_str); 
set(h_Lxt, 'String' ,num2str(lxt)); 
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if Ixt == 0 
Mxt = 1; 

set(findobj(gcf, 'Tag', 'Mxt' ) , 'String', num2str(Mxt)); 
warndlg (' Number of array elements set to 1. ', ' Number of Array 
Elements' , 'help' ); 
end 

Mxt = str2num(get(findobj(gcf,' Tag' , 'Mxt' ), 'String ')) ; 
if Ixt > 0 

if Mxt == 1 

errordlg( 'For the spacing indicated, number of elements must be at 
least 2 . ' , . . . 

'Number of Elements' , 'error' ); 

Ixt = 0; 

set(h_Lxt, 'String' ,num2str (Ixt)); 

end 

end 


case 'Izt' 


h_Lzt = findobj(gcf, 'Tag' , ' Izt' ); 
lzt_str = get(h_Lzt ,' String ') ; 

Izt = getSpacing(lzt_str); 
set(h_Lzt, 'String' ,num2str(Izt)); 
if Izt == 0 
Mzt = 1; 

set(findobj(gcf, 'Tag', 'Mzt'), 'String', num2str(Mzt)); 
warndlg (' Number of array elements set to 1.',... 

' Number of Array Elements ',' warn ') ; 


end 


Mzt = str2num(get(findobj(gcf, 'Tag','Mzt'),'String')); 
if Izt > 0 

if Mzt == 1 

errordlg ('For the spacing indicated, number of elements must be at 
least 2 . ' , . . . 

'Number of Elements' , 'error' ); 

Izt = 0; 

set(h_Lzt, 'String' ,num2str(Izt)); 

end 

end 


case 'txs_xdist' 

h_txs_xdist = findobj(gcf, 'Tag', 'txs_xdist' ); 
txs_xdist_val = get(h_txs_xdist ,' Value ') ; 
ixdist = txs_xdist_val; 
if (ixdist == 2 | ixdist == 4) 

Mxt = str2num(get(findobj(gcf, 'Tag' , 'Mxt' ), 'String' ) ) ; 
if rem(Mxt,2) ~= 0 

set(findobj(gcf, 'Tag' , 'txs_xdist' ), 'Value' ,1); 

errordlg (' Number of subarrays in X-plane must be even. Change 
number of subarrays or select another distribution.', ... 

' Number of Subarrays' , 'error' ); 

end 

end 


if (ixdist == 2) |(ixdist == 3) I(ixdist == 4) 

set(findobj(gcf, 'Tag', 'txs_peddbx '),'Enable','on'); 
set(findobj(gcf, 'Tag', 'txs_nexpx' ), 'Enable' , 'on' ); 
else 

set(findobj(gcf, 'Tag', 'txs_peddbx' ), 'Enable' , 'of f' ); 
set(findobj(gcf, 'Tag', 'txs_nexpx '), 'Enable', 'off'); 

end 
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case 'txs_zdist' 

h_txs_zdist = findobj(gcf, 'Tagtxs_zdist ') ; 
txs_zdist_val = get(h_txs_zdist, 'Value ') ; 
izdist = txs_zdist_val; 
if (izdist == 2 | izdist == 4) 

Mzt = str2num(get(findobj(gcf Tag'MztString ')); 
if rem(Mzt,2) ~= 0 

set(findobj(gcf, 'Tag' , 'txs_zdist' ) , 'Value' ,1); 

errordlg (' Number of subarrays in Z-plane must be even. Change 
number of subarrays or select another distribution.', ... 

' Number of Subarrays' , 'error' ); 

end 

end 

if (izdist == 2) |(izdist == 3) |(izdist == 4) 

set(findobj(gcf, 'Tag', 'txs_peddbz' ), 'Enable' , 'on' ); 
set(findobj(gcf, 'Tag', 'txs_nexpz' ), 'Enable' ,' on' ); 
else 

set(findobj(gcf, 'Tag', 'txs_peddbz' ), 'Enable' , 'off' ); 
set(findobj(gcf, 'Tag', 'txs_nexpz '), 'Enable', 'off'); 

end 

case 'tx_xdist' 

h_tx_xdist = findobj(gcf ,' Tag ',' tx_xdist ') ; 
tx_xdist_val = get(h_tx_xdist,' Value ') ; 
ixdist = tx_xdist_val; 
if (ixdist == 2 | ixdist == 4) 

Nxt = str2num(get(findobj(gcf, 'Tag' , 'Nxt' ), 'String' ) ) ; 
if rem(Nxt,2) ~= 0 

set(findobj(gcf, 'Tag' , 'tx_xdist' ),' Value' ,1); 

errordlg (' Number of elements in X-plane must be even. Change number 
of elements or select another distribution.', ... 

' Number of Subrray Elements ',' error ') ; 

end 

end 

if (ixdist == 2) |(ixdist == 3) |(ixdist == 4) 

set(findobj(gcf, 'Tag', 'tx_peddbx' ), 'Enable' ,' on' ); 
set(findobj(gcf, 'Tag', 'tx_nexpx' ),' Enable' , 'on' ); 
else 

set(findobj(gcf, 'Tag', 'tx_peddbx '), 'Enable', 'off'); 
set(findobj(gcf, 'Tag', 'tx_nexpx' ),' Enable' , 'off' ); 

end 

case 'tx_zdist' 

h_tx_zdist = findobj(gcf ,' Tag ',' tx_zdist ') ; 
tx_zdist_val = get(h_tx_zdist,' Value ') ; 
izdist = tx_zdist_val; 
if (izdist == 2 | izdist == 4) 

Nzt = str2num(get(findobj(gcf ,' Tag ',' Nzt '),' String ')) ; 
if rem(Nzt,2) ~= 0 

set(findobj(gcf, 'Tag' , 'tx_zdist' ),' Value' ,1); 

errordlg (' Number of elements in Z-plane must be even. Change number 
of elements or select another distribution.', ... 

' Number of Subrray Elements ',' error ') ; 

end 

end 

if (izdist == 2) |(izdist == 3) |(izdist == 4) 

set(findobj(gcf, 'Tag', 'tx_peddbz' ), 'Enable' ,' on' ); 
set(findobj(gcf, 'Tag', 'tx_nexpz' ),' Enable' , 'on' ); 
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else 

set(findobj(gcf, 'Tag', 'tx_peddbz '), 'Enable', 'off'); 
set(findobj(gcf, 'Tag' , 'tx_nexpz' ), 'Enable' , 'off' ); 

end 


%Case for Receive Array parameters 
case 'Nxr' 


zero 


h_Nxr = findobj(gcf, 'Tag ',' Nxr ') ; 
nxr_str = get(h_Nxr, 'String'); 

Nxr = getNEL(nxr_str); 
set(h_Nxr, 'String' ,num2str(Nxr)); 
if Nxr == 1 
dxr = 0; 

set(findobj(gcf, 'Tag','dxr'), 'String', num2str(dxr)); 
warndlg('For a single element, array spacing is set to zero.',... 
' Array Spacing ',' help '); 
elseif Nxr > 1 

dxr = str2num(get(findobj(gcf ,' Tag ',' dxr '),' String ')); 
if dxr == 0 

set(findobj(gcf, 'Tag' , 'dxr' ), 'String' ,num2str(1)); 
warndlg('For the number of elements, array spacing cannot be 

t • • • 

' Array Spacing ',' help ') ; 

end 


end 


case 'Nzr' 

h_Nzr = findobj(gcf,' Tag' ,' Nzr' ); 
nzr_str = get(h_Nzr,' String ') ; 

Nzr = getNEL(nzr_str); 
set(h_Nzr,' String' ,num2str(Nzr)); 
if Nzr == 1 
dzr = 0; 

set(findobj(gcf, 'Tag','dzr'),'String', num2str(dzr)); 
warndlgt'For a single element, array spacing is set zero.',... 

' Array Spacing ',' help '); 
elseif Nzr > 1 

dzr = str2num(get(findobj(gcf ,' Tag ',' dzr '),' String ')) ; 
if dzr == 0 

set(findobj(gcf, 'Tag', 'dzr'), 'String', num2str(1)); 
warndlg('For the number of elements, array spacing cannot be 

zero .',... 


end 


Array Spacing ','help'); 


end 


case 'dxr' 

h_Dxr = findobj(gcf,' Tag ',' dxr ') ; 
dxr_str = get(h_Dxr,' String ') ; 
dxr = getSpacing(dxr_str); 
set(h_Dxr,' String' ,num2str(dxr)); 
if dxr == 0 
Nxr = 1; 

set(findobj(gcf, 'Tag', 'Nxr' ), 'String' ,num2str(Nxr)); 
warndlg (' Number of array elements set to 1. ', ' Number of Array 
Elements' , 'help' ); 
end 


Nxr = str2num(get(findobj(gcf, 'Tag' ,' Nxr' ),' String' )); 
if dxr > 0 

if Nxr == 1 
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errordlg ( 'For the spacing indicated, number of elements must be at 
least 2 . ' , . . . 

'Number of Elements' , 'error' ); 
dxr = 0; 

set(h_Dxr, ' String' , num2str(dxr)); 

end 

end 


case 'dzr' 


h_Dzr = findobj(gcf, 'Tag' ,' dzr' ); 
dzr_str = get(h_Dzr, 'String' ); 
dzr = getSpacing(dzr_str); 
set(h_Dzr, 'String' ,num2str(dzr)); 
if dzr == 0 
Nzr = 1; 

set(findobj(gcf, 'Tag','Nzr'),'String', num2str(Nzr)); 
warndlg( 'Number of array elements set to 1.' , ... 

' Number of Array Elements ',' warn ') ; 


end 


Nzr = str2num(get(findobj(gcf,' Tag' ,' Nzr' ),' String' ) ) ; 
if dzr > 0 

if Nzr == 1 

errordlg ('For the spacing indicated, number of elements must be at 
least 2.' ,. .. 

'Number of Elements' , 'error' ); 
dzr = 0; 

set(h_Dzr,' String' ,num2str(dzr)); 

end 

end 


case 'Mxr' 


zero 


h_Mxr = findobj(gcf ,' Tag ',' Mxr ') ; 
mxr_str = get(h_Mxr, 'String'); 

Mxr = getNEL(mxr_str); 
set(h_Mxr, 'String' ,num2str(Mxr)); 
if Mxr == 1 
Ixr = 0; 

set(findobj(gcf, 'Tag', 'lxr' ), 'String', num2str(Ixr)); 
warndlg ('For a single element, array spacing is set to zero.',... 
' Array Spacing ',' help '); 
elseif Mxr > 1 

lxr = str2num(get(findobj(gcf ,' Tag ',' lxr '),' String ')); 
if lxr == 0 

set(findobj(gcf, 'Tag', 'lxr'), 'String', num2str(1)); 
warndlg('For the number of elements, array spacing cannot be 

r . . . 

' Array Spacing ',' help ') ; 

end 


end 


case 'Mzr' 

h_Mzr = findobj(gcf, 'Tag' , 'Mzr' ); 
mzr_str = get(h_Mzr, 'String'); 

Mzr = getNEL(mzr_str); 
set(h_Mzr, 'String' ,num2str(Mzr)); 
if Mzr == 1 
1 z r = 0; 

set(findobj(gcf, 'Tag', '1zr' ), 'String', num2str(Izr)); 
warndlg ('For a single element, array spacing is set zero.',... 
' Array Spacing ',' help '); 
elseif Mzr > 1 
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zero 


lzr = str2num(get(findobj(gcf, 'TagIzrString ')) ; 
if lzr == 0 

set(findobj(gcf, 'Tag 1 ,'lzr'),'String', num2str(1)); 
warndlgf 'For the number of elements, array spacing cannot be 

/ • • • 

' Array Spacing ',' help ') ; 

end 

end 


case 'Ixr' 

h_Lxr = findobj(gcf,' Tag' , 'Ixr' ); 
lxr_str = get(h_Lxr,' String ') ; 

Ixr = getSpacing(lxr_str); 
set(h_Lxr,' String' ,num2str(Ixr)); 
if Ixr == 0 
Mxr = 1; 

set(findobj(gcf, 'Tag', 'Mxr' ), 'String' ,num2str(Mxr)); 
warndlg (' Number of array elements set to 1. ', ' Number of Array 
Elements' , 'help' ); 
end 

Mxr = str2num(get(findobj(gcf,' Tag' , 'Mxr' ),' String' )); 
if Ixr > 0 

if Mxr == 1 

errordlg( 'For the spacing indicated, number of elements must be at 
least 2.' ,... 

'Number of Elements' , 'error' ); 

Ixr = 0; 

set(h_Lxr, 'String' ,num2str(Ixr)); 

end 

end 


case 'lzr' 


h_Lzr = findobj(gcf, 'Tag' , ' lzr ' ); 
lzr_str = get(h_Lzr, 'String ') ; 
lzr = getSpacing(lzr_str); 
set(h_Lzr, 'String' ,num2str(lzr)); 
if lzr == 0 
Mzr = 1; 

set(findobj(gcf, 'Tag','Mzr'), 'String', num2str(Mzr)); 
warndlg( 'Number of array elements set to 1 .' , ... 

' Number of Array Elements ',' warn ') ; 


end 


Mzr = str2num(get(findobj(gcf,' Tag' , 'Mzr' ),' String' )); 
if lzr > 0 

if Mzr == 1 

errordlg ('For the spacing indicated, number of elements must be at 
least 2 . ' , . . . 

'Number of Elements' , 'error' ); 
lzr = 0; 

set(h_Lzr,' String' ,num2str(lzr)); 

end 

end 


case 'rxs_xdist' 

h_rxs_xdist = findobj(gcf, 'Tag', 'rxs_xdist' ); 
rxs_xdist_val = get(h_rxs_xdist ,' Value ') ; 
ixdist = rxs_xdist_val; 
if (ixdist == 2 | ixdist == 4) 

Mxr = str2num(get(findobj(gcf, 'Tag' , 'Mxr' ), 'String' )); 
if rem(Mxr,2) ~= 0 
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set(findobj(gcf, 'Tag' , 'rxs_xdist' ), 'Value' ,1); 

errordlg (' Number of subarrays in X-plane must be even. Change 
number of subarrays or select another distribution.', ... 

' Number of Subarrays' , 'error' ); 

end 

end 

if (ixdist == 2) |(ixdist == 3) |(ixdist == 4) 

set(findobj(gcf, 'Tag', 'rxs_peddbx' ), 'Enable' , 'on' ); 
set(findobj(gcf, 'Tag', 'rxs_nexpx' ), 'Enable' ,' on' ); 
else 

set(findobj(gcf, 'Tag', 'rxs_peddbx' ), 'Enable' , 'off' ); 
set(findobj(gcf, 'Tag', 'rxs_nexpx '), 'Enable','off'); 

end 

case 'rxs_zdist' 

h_rxs_zdist = findobj(gcf, 'Tag' , 'rxs_zdist' ); 
rxs_zdist_val = get(h_rxs_zdist, 'Value ') ; 
izdist = rxs_zdist_val; 
if (izdist == 2 | izdist == 4) 

Mzr = str2num(get(findobj(gcf, 1 Tag' , ' Mzr' ), ' String' )); 
if rem(Mzr,2) ~= 0 

set(findobj(gcf, 'Tag' , 'rxs_zdist' ) , 'Value' ,1); 

errordlg (' Number of subarrays in Z-plane must be even. Change 
number of subarrays or select another distribution.', ... 

' Number of Subarrays' , 'error' ); 

end 

end 


if (izdist == 2) 
set(findobj(gcf, 
set (findobj(gcf, 
else 

set(findobj(gcf, 
set(findobj(gcf, 
end 


(izdist == 3) |(izdist == 4) 

'Tag' , 'rxs_peddbz' ), 'Enable' , ' on ' ); 

'Tag' , 'rxs_nexpz' ), 'Enable' , 'on' ); 

'Tag' , 'rxs_peddbz' ), 'Enable' , 'off' ); 
'Tag' , 'rxs_nexpz' ), 'Enable' , 'of f' ); 


case 'rx_xdist' 

h_rx_xdist = findobj(gcf, 'Tag',' rx_xdist' ); 
rx_xdist_val = get(h_rx_xdist ,' Value ') ; 
ixdist = rx_xdist_val; 
if (ixdist == 2 | ixdist == 4) 

Nxr = str2num(get(findobj(gcf, 'Tag', 'Nxr' ), 'String' ) ) ; 
if rem(Nxr,2) ~= 0 

set(findobj(gcf, 'Tag' , 'rx_xdist' ),' Value' , 1) ; 

errordlg (' Number of elements in X-plane must be even. Change number 
of elements or select another distribution.', ... 

' Number of Subrray Elements ',' error ') ; 

end 

end 

if (ixdist == 2) |(ixdist == 3) I(ixdist == 4) 

set(findobj(gcf, 'Tag' ,' rx_peddbx' ), 'Enable' , 'on' ); 
set(findobj(gcf, 'Tag', 'rx_nexpx' ), 'Enable' , 'on' ); 
else 

set(findobj(gcf, 'Tag' ,' rx_peddbx '), 'Enable','off'); 
set(findobj(gcf, 'Tag' ,' rx_nexpx' ), 'Enable' , 'off' ); 

end 


case 'rx_zdist' 

h_rx_zdist = findobj(gcf, 'Tag','rx_zdist'); 
rx_zdist_val = get(h_rx_zdist, 'Value ') ; 
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izdist = rx_zdist_val; 
if (izdist == 2 | izdist == 4) 

Nzr = str2num(get(findobj(gcf, 'Tag', 'Nzr' ), 'String' )); 
if rem(Nzr,2) ~= 0 

set(findobj(gcf, 'Tag' , 'rx_zdist' ),' Value' ,1); 

errordlg (' Number of elements in Z-plane must be even. Change number 
of elements or select another distribution.', ... 

' Number of Subrray Elements','error'); 

end 

end 

if (izdist == 2) |(izdist == 3) I(izdist == 4) 

set(findobj(gcf, 'Tag' , 'rx_peddbz' ),' Enable' , 'on' ); 
set(findobj(gcf, 'Tag' ,' rx_nexpz '),'Enable','on'); 
else 

set(findobj(gcf, 'Tag' , 'rx_peddbz' ),' Enable' , 'off' ); 
set(findobj(gcf, 'Tag', 'rx_nexpz' ), 'Enable' , 'off' ); 

end 

%Case Calculate and plot two-way pattern 
case 'Calculate' 

rad=pi/180; 

%Get parameters 

h_eltype = get(findobj(gcf ,' Tag ',' eltype '),' Value ') ; 


if 

h_eltype 

== 1, 

diptype = 

' h' ; 

: end 

if 

h_eltype 

== 2, 

diptype = 

's'; 

: end 

if 

h_eltype 

== 3, 

diptype = 

' i' J 

: end 


h_dipdir = get(findobj(gcf ,' Tagdipdir '), 'Value ') ; 
if h_dipdir == 1, dipdir = 'x'; end 
if h_dipdir == 2, dipdir = 'z'; end 
f = (str2num(get(findobj(gcf,'Tag','frequency'),'String')))*le6; 
h = str2num(get(findobj(gcf ,' Tag ',' height '),' String ')) ; 
thetas = str2num(get(findobj(gcf, 'Tag',' thetas' ), 'String ')) ; 
phis = str2num(get(findobj(gcf ,' Tag ',' phis '),' String ')) ; 
tstart = str2num(get(findobj(gcf, 'Tag',' tstart ') , 'String')); 
tstop = str2num(get(findobj(gcf ,' Tag ',' tstop '),' String ')) ; 
pstart = str2num(get(findobj(gcf, 'Tag',' pstart' ), 'String')); 


pstop 

= str2num(get(findobj(gcf,'Tag','pstop'),'String 

delt 

= str2num(get(findobj(gcf 

, 'Tag' 

, 1 delt' 

, 'String 

) 

delp 

= str2num(get(findobj(gcf 

, 'Tag' 

, 'delp' 

, 'String 

) 

Nxt = 

str2num(get 

(findobj(gcf. 

'Tag' , 

’Nxt' ), 

'String' ) 

r 

Nzt = 

str2num(get 

(findobj(gcf. 

’Tag' , 

1 Nzt' ), 

'String' ) 

f 

dxt = 

str2num(get 

(findobj(gcf. 

’Tag' , 

'dxt' ) , 

'String' ) 

r 

dzt = 

str2num(get 

(findobj(gcf. 

’Tag' , 

'dzt' ) , 

'String' ) 

f 

Mxt = 

str2num(get 

(findobj(gcf. 

’Tag' , 

'Mxt' ), 

'String' ) 

r 

Mzt = 

str2num(get 

(findobj(gcf. 

'Tag' , 

'Mzt' ), 

'String' ) 

r 

lxt = 

str2num(get 

(findobj(gcf. 

’Tag' , 

'lxt' ) , 

'String' ) 

r 

lzt = 

str2num(get 

(findobj(gcf. 

’Tag' , 

'lzt' ) , 

'String' ) 

r 

Nxr = 

str2num(get 

(findobj(gcf. 

’Tag' , 

’Nxr ' ), 

'String' ) 

r 

Nzr = 

str2num(get 

(findobj(gcf. 

'Tag' , 

’Nzr' ), 

'String' ) 

r 

dxr = 

str2num(get 

(findobj(gcf. 

’Tag' , 

'dxr' ), 

'String' ) 

r 

dzr = 

str2num(get 

(findobj(gcf. 

’Tag' , 

’dzr 1 ), 

'String' ) 

r 

Mxr = 

str2num(get 

(findobj(gcf. 

'Tag' , 

'Mxr' ), 

'String' ) 

r 

Mzr = 

str2num(get 

(findobj(gcf. 

’Tag' , 

'Mzr' ), 

'String' ) 

r 

lxr = 

str2num(get 

(findobj(gcf. 

’Tag' , 

'lxr ' ), 

'String' ) 

f 

lzr = 

str2num(get 

(findobj(gcf. 

’Tag' , 

'lzr'), 

'String' ) 

r 


txs_xdist = get(findobj(gcf ,' Tag ',' txs_xdist '),' Value ') ; 
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txs_peddbxval = get(findobj(gcf, 'Tagtxs_peddbx '), 'Value' ) 


if txs_peddbxval==l, 
if txs_peddbxval==2, 
if txs peddbxval==3, 
if txs peddbxval==4, 
if txs_peddbxval==5, 
if txs_peddbxval==6, 
if txs peddbxval==7, 
if txs peddbxval==8, 


txs_peddbx=15; end 
txs_peddbx=20; end 
txs_peddbx=25; end 
txs_peddbx=30; end 
txs_peddbx=35; end 
txs_peddbx=40; end 
txs_peddbx=45; end 
txs_peddbx=50; end 


txs_nexpx = get(findobj(gcf Tagtxs_nexpxValue ') ; 

txs_zdist = get(findobj(gcf, 'Tagtxs_zdistValue ’) ; 
txs_peddbzval = get(findobj(gcf, 'Tagtxs_peddbz' ), 'Value' ) 
if txs_peddbzval==l, txs_peddbz=15; end 

txs_peddbz=20; end 
txs_peddbz=25; end 
txs_peddbz=30; end 
txs_peddbz=35; end 
txs_peddbz=40; end 
txs_peddbz=45; end 
txs_peddbz=50; end 


if txs_peddbzval==2, 
if txs peddbzval==3, 
if txs peddbzval==4, 
if txs_peddbzval==5, 
if txs_peddbzval==6, 
if txs peddbzval==7, 
if txs_peddbzval==8, 


txs_nexpz = get(findobj(gcf, 'Tag ',' txs_nexpz '),' Value ') ; 

tx_xdist = get(findobj(gcf,' Tag ',' tx_xdist '), 'Value ') ; 
tx_peddbxval = get(findobj(gcf, 'Tagtx_peddbx '), 'Value'); 


if 

tx_peddbxval= 

=1, 

tx_peddbx=15; 

end 

if 

tx_peddbxval= 

=2, 

tx_peddbx=20; 

end 

if 

tx_peddbxval= 

= 3, 

tx_peddbx=25; 

end 

if 

tx_peddbxval= 

= 4, 

tx_peddbx=30 ; 

end 

if 

tx_peddbxval= 

= 5, 

tx_peddbx=35 ; 

end 

if 

tx_peddbxval= 

= 6, 

tx_peddbx=40; 

end 

if 

tx_peddbxval= 

= 7, 

tx_peddbx=45; 

end 

if 

tx_peddbxval= 

= 8, 

tx_peddbx=50; 

end 


tx_nexpx = get(findobj(gcf Tag ',' tx_nexpx '),' Value ') ; 

tx_zdist = get(findobj(gcf,' Tag ',' tx_zdist '), 'Value ') ; 
tx_peddbzval = get(findobj(gcf ,' Tagtx_peddbz '), 'Value') ; 
if tx_peddbzval==l, tx peddbz=15; end 
tx_peddbz=20; end 
tx_peddbz=25; end 
tx_peddbz=30; end 
tx_peddbz=35; end 
tx_peddbz=40; end 
tx_peddbz=45; end 
tx_peddbz=50; end 


if tx_peddbzval==2, 
if tx_peddbzval==3, 
if tx_peddbzval==4, 
if tx_peddbzval==5, 
if tx_peddbzval==6, 
if tx_peddbzval==7, 
if tx_peddbzval==8, 


tx_nexpz = get(findobj(gcf,' Tag ',' tx_nexpz '),' Value ') ; 


txs_ampx = getamplitudes(Mxt,txs_xdist,txs_peddbx,txs_nexpx) 
txs_ampz = getamplitudes(Mzt,txs_zdist,txs_peddbz,txs_nexpz) 
tx_ampx = getamplitudes(Nxt,tx_xdist,tx_peddbx,tx_nexpx); 
tx_ampz = getamplitudes(Nzt,tx_zdist,tx_peddbz,tx_nexpz); 


rxs_xdist = get(findobj(gcf ,' Tag ',' rxs_xdist '),' Value ') ; 
rxs_peddbxval = get(findobj(gcf,' Tag ',' rxs_peddbx '), 'Value' ) 


if 

rxs peddbxval= 

= 1, 

rxs_peddbx=15; 

end 

if 

rxs_peddbxval= 

=2, 

rxs peddbx=20; 

end 

if 

rxs_peddbxval= 

=3, 

rxs peddbx=25; 

end 

if 

rxs_peddbxval= 

= 4, 

rxs_peddbx=30 ; 

end 

if 

rxs_peddbxval= 

=5, 

rxs_peddbx=35 ; 

end 

if 

rxs_peddbxval= 

= 6, 

rxs_peddbx=40; 

end 

if 

rxs_peddbxval= 

=7, 

rxs_peddbx=45; 

end 

if 

rxs peddbxval= 

= 8, 

rxs_peddbx=50; 

end 
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rxs_nexpx = get(findobj(gcf, 'Tagrxs_nexpxValue ') ; 

rxs_zdist = get(findobj(gcf Tagrxs_zdistValue ') ; 
rxs_peddbzval = get(findobj(gcf,' Tag ',' rxs_peddbz '), 'Value ') ; 


if 

rxs_peddbzval= 

= 1, 

rxs_peddbz=15; 

end 

if 

rxs_peddbzval= 

=2, 

rxs peddbz=20; 

end 

if 

rxs_peddbzval= 

=3, 

rxs peddbz=25; 

end 

if 

rxs_peddbzval= 

= 4, 

rxs_peddbz=30; 

end 

if 

rxs_peddbzval= 

=5, 

rxs_peddbz=35 ; 

end 

if 

rxs_peddbzval= 

= 6, 

rxs_peddbz=40; 

end 

if 

rxs_peddbzval= 

=7, 

rxs_peddbz=45; 

end 

if 

rxs_peddbzval= 

= 8, 

rxs_peddbz=50; 

end 


rxs_nexpz = get(findobj(gcf Tag 1 rxs_nexpzValue ') ; 


rx_xdist = get(findobj(gcf Tagrx_xdist '), 'Value ') ; 
rx_peddbxval = get(findobj(gcf ,' Tagrx_peddbx '), 'Value '); 
if rx_peddbxval==l, rx peddbx=15; end 
rx_peddbx=20; end 
rx_peddbx=25; end 
rx_peddbx=30; end 
rx_peddbx=35; end 
rx_peddbx=40; end 
rx_peddbx=45; end 
rx_peddbx=50; end 


if 

if 

if 

if 

if 

if 

if 


rx_peddbxval==2, 
rx_peddbxval==3, 
rx_peddbxval==4, 
rx_peddbxval==5, 
rx_peddbxval==6, 
rx_peddbxval==7, 
rx_peddbxval==8, 


rx_nexpx = get(findobj(gcf,' Tag ',' rx_nexpx '),' Value ') ; 


rx_zdist = get(findobj(gcf, 'Tag ',' rx_zdist '), 'Value ') ; 


eddbzval = get(findobj(gcf ,' Tag' , 

'rx_peddbz 

if 

rx_peddbzval= 

=1, 

rx_peddbz=15; 

end 

if 

rx_peddbzval= 

=2, 

rx_peddbz=20; 

end 

if 

rx_peddbzval= 

= 3, 

rx_peddbz=25; 

end 

if 

rx_peddbzval= 

= 4, 

rx_peddbz=30 ; 

end 

if 

rx_peddbzval= 

=5, 

rx_peddbz=35 ; 

end 

if 

rx_peddbzval= 

= 6, 

rx_peddbz=40; 

end 

if 

rx_peddbzval= 

= 7, 

rx_peddbz=45; 

end 

if 

rx_peddbzval= 

= 8, 

rx_peddbz=50; 

end 


rx_nexpz = get(findobj(gcf ,' Tagrx_nexpz '), 'Value ') ; 


rxs_ampx = getamplitudes(Mxr,rxs_xdist,rxs_peddbx,rxs_nexpx); 
rxs_ampz = getamplitudes(Mzr,rxs_zdist,rxs_peddbz,rxs_nexpz); 
rx_ampx = getamplitudes(Nxr,rx_xdist,rx_peddbx,rx_nexpx); 
rx_ampz = getamplitudes(Nzr,rx_zdist,rx_peddbz,rx_nexpz); 


% Compute CAF of Transmit DSA 

if diptype == 'h' %compute for half wave dipoles 

[TxSAFtheta,TxSAFphi,TxDAFtheta,TxDAFphi,Ut,Wt]=caf_hdip2( 't' ,f,dxt,Nxt,dzt,Nzt 
,Ixt,Mxt,Izt,Mzt,h,dipdir,thetas,phis,tstart,tstop,pstart,pstop,delt,delp,txs_a 
mpx,txs_ampz,tx_ampx,tx_ampz); 
end 


if diptype == 's' %compute for short dipoles 

[TxSAFtheta,TxSAFphi,TxDAFtheta,TxDAFphi,Ut,Wt]=caf_sdip2( 't' ,f,dxt,Nxt,dzt,Nzt 
,Ixt,Mxt,Izt,Mzt,h,dipdir,thetas,phis,tstart,tstop,pstart,pstop,delt,delp, txs_a 
mpx,txs_ampz,tx_ampx,tx_ampz); 
end 


if diptype == 'i' %compute for isotropic elements 
[TxSAFtheta,TxSAFphi,TxDAFtheta,TxDAFphi,Ut,Wt]=caf_iso2( 't' ,f,dxt,Nxt,dzt,Nzt, 
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lxt,Mxt, lzt,Mzt,h,thetas,phis,tstart,tstop,pstart,pstop,delt,delp,txs_ampx,txs_ 
arapz,tx_ampx,tx_ampz); 
end 

% Normalize transmit DSA pattern 

MAXtx=max([max(max(abs(TxDAFtheta))),max(max(abs(TxDAFphi)))]); 
TxDAFt_norm=abs(TxDAFtheta)/MAXtx; %Find magnitude and normalize 
TxDAFp_norm=abs(TxDAFphi)/MAXtx; 

TxDAFt_db=20*logl0(TxDAFt_norm); %Compute in dB 
TxDAFp_db=20*logl0(TxDAFp_norm); 

% Compute CAF of Receive DSA 

if diptype == 'h' %Compute for half wave dipoles 

[RxSAFtheta,RxSAFphi,RxDAFtheta,RxDAFphi,Ur,Wr]=caf_hdip2( 1 r' ,f,dxr,Nxr,dzr,Nzr 
,lxr,Mxr,lzr,Mzr,h,dipdir,thetas,phis,tstart,tstop,pstart,pstop,delt,delp,rxs_a 
mpx,rxs_ampz,rx_ampx,rx_ampz); 
end 

if diptype == 's' %Compute for short dipoles 

[RxSAFtheta,RxSAFphi,RxDAFtheta,RxDAFphi,Ur,Wr]=caf_sdip2( 'r' ,f,dxr,Nxr,dzr,Nzr 
, lxr,Mxr, lzr,Mzr,h, dipdir,thetas,phis,tstart,tstop,pstart,pstop,delt,delp,rxs_a 
mpx,rxs_ampz,rx_ampx,rx_ampz); 
end 

if diptype == 'i' %Compute for isotropic elements 

[RxSAFtheta,RxSAFphi,RxDAFtheta,RxDAFphi,Ur,Wr]=caf_iso2('r ’ ,f,dxr,Nxr,dzr,Nzr, 
lxr,Mxr,lzr,Mzr,h,thetas,phis,tstart,tstop,pstart,pstop,delt,delp,rxs_ampx,rxs_ 
ampz,rx_ampx,rx_ampz); 
end 

% Normalize receive DSA pattern 

MAXrx=max([max(max(abs(RxDAFtheta))),max(max(abs(RxDAFphi)))]); 
RxDAFt_norm=abs(RxDAFtheta)/MAXrx; %Find magnitude and normalize 
RxDAFp_norm=abs(RxDAFphi)/MAXrx; 

RxDAFt_db=20*logl0(RxDAFt_norm); %Compute in dB 
RxDAFp_db=20*logl0(RxDAFp_norm); 

% Compute Two-way Pattern of Transmit and Receive DSAs 
Two_way_DAFt=TxDAFtheta.*RxDAFtheta; 

Two_way_DAFp=TxDAFphi.*RxDAFphi; 

MAXtw=max([max(max(abs(Two_way_DAFt))),max(max(abs(Two_way_DAFp)))]); 
Two_way_DAFt_norm=abs(Two_way_DAFt)/MAXtw; %Normalized theta pattern 
Two_way_DAFp_norm=abs(Two_way_DAFp)/MAXtw; %Normalized phi pattern 

Two_way_DAFt_db=20*logl0(Two_way_DAFt_norm) ; %Compute in dB 
Two_way_DAFp_db=20*logl0(Two_way_DAFp_norm); 

%Export dsa configuration to file dsaconfig.m 
save dsaconfig f diptype dipdir h delt delp thetas phis ... 
tstart tstop pstart pstop ... 

Nxt Nzt dxt dzt Mxt Mzt lxt lzt . . . 

Nxr Nzr dxr dzr Mxr Mzr lxr lzr ... 
txs_xdist txs_peddbx txs_nexpx ... 
txs_zdist txs_peddbz txs_nexpz ... 
tx_xdist tx_peddbx tx_nexpx ... 
tx_zdist tx_peddbz tx_nexpz ... 
rxs_xdist rxs_peddbx rxs_nexpx ... 
rxs_zdist rxs_peddbz rxs_nexpz ... 
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rx_xdist rx_peddbx rx_nexpx ... 
rx_zdist rx_peddbz rx_nexpz ... 

%Export dsa pattern data to file dsapattern.m 

save dsapattern Ur Wr TxDAFt_db TxDAFp_db RxDAFt_db RxDAFp_db ... 
Two_way_DAFt_db Two_way_DAFp_db 

%Set plot dynamic range to pmin 
pmin = -80; %Set the pmin value (dB) 
ip = floor((pstop-pstart)/delp)+1; 
it = floor((tstop-tstart)/delt)+1; 
for il = 1:ip 

for i2 = 1:it 

theta(il,i2) = tstart + (i2-l)*delt; 
phi(il,i2) = pstart + (il-l)*delp; 

end 

end 

if ip == 1 %Plot phi cut 
figure (1),elf 
subplot(211) 

plot(theta,TxDAFt_db),grid %Plot transmit DSA theta component 
axis([tstart,tstop,pmin,0]) 

title (['Transmit DSA Pattern, \phi= ', num2str(pstart ),' A o cut']) 
xlabel (' Pattern Angle, \theta (deg)') 
ylabel (' Normalized F_\thetaI (dB) ') 
subplot(212) 

plot(theta,TxDAFp_db),grid %Plot transmit DSA phi component 
axis ( [tstart,tstop,pmin,0]) 

title (['Transmit DSA Pattern, \phi= ', num2str(pstart ),' A o cut']) 
xlabel (' Pattern Angle, \theta (deg)') 
ylabel( 'Normalized F_\phi| (dB) ') 


figure(2),elf 
subplot(211) 

plot(theta,RxDAFt_db),grid %Plot receive DSA theta component 
axis ( [tstart,tstop,pmin,0]) 

title ([ Receive DSA Pattern, \phi= ',num2str(pstart) ,' A o cut']) 
xlabel (' Pattern Angle, \theta (deg)') 
ylabel(' Normalized |F_\theta| (dB) ') 
subplot(212) 

plot(theta,RxDAFp_db),grid %Plot receive DSA phi component 
axis ( [tstart,tstop,pmin,0]) 

title ( [ Receive DSA Pattern, \phi= ',num2str(pstart) ,' A o cut']) 
xlabel (' Pattern Angle, \theta (deg)') 
ylabel(' Normalized F_\phi| (dB) ') 

figure(3),elf 
subplot(211) 

plot(theta,Two_way_DAFt_db),grid %Plot two-way theta component 
axis([tstart,tstop,pmin,0]) 

title ([ Two-way Pattern, \phi= ',num2str(pstart) ,' A o cut']) 
xlabel (' Pattern Angle, \theta (deg)') 
ylabel(' Normalized F_\theta| (dB) ') 
subplot(212) 

plot(theta,Two_way_DAFp_db),grid %Plot two-way phi component 
axis([tstart,tstop,pmin,0]) 

title ( [ Two-way Pattern IF_n_o_r_m(\phi) |, 

\phi=' ,num2str(pstart), ' A o cut' ]) 

xlabel (' Pattern Angle, \theta (deg)') 
ylabel(' Normalized F_\phi| (dB) ') 
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end 


if it == 1 %Plot theta cut 
figure (1),elf 
subplot(211) 

plot(phi,TxDAFt_db),grid %Plot transmit DSA theta component 
axis([pstart,pstop,pmin,0]) 

title (['Transmit DSA Pattern, \theta= ', num2str(tstart ),' A o cut']) 
xlabel (' Pattern Angle, \phi (deg) 1 ) 
ylabel (' Normalized |F_\theta| (dB) ') 
subplot(212) 

plot(phi,TxDAFp_db),grid %Plot transmit DSA phi component 
axis ( [pstart,pstop,pmin,0]) 

title ([' Transmit DSA Pattern, \theta= ', num2str(tstart ),' A o cut']) 
xlabel (' Pattern Angle, \phi (deg) ') 
ylabel (' Normalized IF_\phi| (dB) ') 

figure(2),elf 
subplot(211) 

plot(phi,RxDAFt_db),grid %Plot receive DSA theta component 
axis ( [pstart,pstop,pmin,0]) 

title ([' Receive DSA Pattern, \theta= ', num2str(tstart) ,' A o cut']) 
xlabel (' Pattern Angle, \phi (deg) ') 
ylabel( 'Normalized |F_\theta| (dB) ') 
subplot(212) 

plot(phi,RxDAFp_db),grid %Plot receive DSA phi component 
axis ( [pstart,pstop,pmin,0]) 

title ([' Receive DSA Pattern, \theta= ', num2str(tstart) ,' A o cut ]) 
xlabel (' Pattern Angle, \phi (deg)') 
ylabel (' Normalized F_\phi| (dB) ') 

figure(3),elf 
subplot(211) 

plot(phi,Two_way_DAFt_db),grid %Plot two-way theta component 
axis([pstart,pstop,pmin,0]) 

title ([' Two-way Pattern, \theta= ',num2str(tstart) ,' A o cut ]) 
xlabel (' Pattern Angle, \phi (deg)') 
ylabel (' Normalized F_\theta| (dB) ') 
subplot(212) 

plot(phi,Two_way_DAFp_db),grid %Plot two-way phi component 
axis([pstart,pstop,pmin,0]) 

title ([' Two-way Pattern, \theta= ', num2str(tstart) ,' A o cut ]) 
xlabel (' Pattern Angle, \phi (deg) ') 
ylabel( 'Normalized IF_\phi| (dB) ') 

end 


if ip>l & it>l 


for il = 
for 


1:ip %Set 
i2 = 1: it 
if TxDAFt 
TxDAFt 
if TxDAFp_i 
TxDAFp_i 
if RxDAFt 
RxDAFt 
if RxDAFp_i 
RxDAFp_i 
if Two_way. 

Two_way. 
if Two_way. 
Two_way. 


mesh plot dynamic range to pmin 


_db (il, i2) < pmin, . . . 
_db(il,i2) = pmin; end 
db (il,i2) < pmin, ... 
db(il,i2) = pmin; end 
db (il,i2) < pmin, .. . 
_db(il,i2) = pmin; end 
.db(il,i2) < pmin, ... 
db(il,i2) = pmin; end 
DAFt_db(il,i2) < pmin,... 

_DAFt_db(il,i2) = pmin; end 
■_DAFp_db(il,i2) < pmin, 

■_DAFp_db(il,i2) = pmin; 
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end 


end 


end 

figure (1),elf %Plot theta component of transmit array 
meshc(Ut,Wt,TxDAFt_db),grid,axis([-1 1-11 pmin,0]),grid 
axis square 

xlabel('U = sin\theta*cos\phi' ) 
ylabel('W = cosXtheta') 
zlabel( 'Normalized F_\theta| (dB) ') 
view (45, 45) 

title( 1 Transmit DSA Pattern') 

figure(2),elf %Plot phi component of transmit array 
mesh(Ut,Wt,TxDAFp_db),grid,axis([-1 1-11 pmin,0]),grid 
axis square 

xlabel('U = sin\theta*cos\phi' ) 
ylabel('W = cosXtheta') 
zlabel( 'Normalized IF_\phi| (dB) ') 
view (45, 45) 

title(' Transmit DSA Pattern') 

figure(3),elf %Plot theta component of receive array 
meshc(Ur,Wr,RxDAFt_db),grid,axis([-1 1-11 pmin,0]),grid 
axis square 

xlabel('U = sin\theta*cos\phi' ) 
ylabel('W = cosXtheta') 
zlabel( 'Normalized F_\thetaI (dB) ') 
view (45, 45) 

title (' Receive DSA Pattern') 

figure(4),elf %Plot phi component of receive array 
meshc(Ur,Wr,RxDAFp_db),grid,axis([-1 1-11 pmin,0]),grid 
axis square 

xlabel('U = sin\theta*cos\phi' ) 
ylabel('W = cosXtheta') 
zlabel( 'Normalized IF_\phi| (dB) ') 
view (45, 45) 

title (' Receive DSA Pattern') 

figure(5),elf %Plot theta component of two-way pattern 

meshc(Ur,Wr,Two_way_DAFt_db),grid,axis([-1 1-11 pmin,0]),grid 

axis square 

xlabel('U = sin\theta*cos\phi' ) 
ylabel('W = cosXtheta') 
zlabel (' Normalized |F_\theta| (dB) ') 
view (45, 45) 

title (' Two-way Pattern') 

figure(6),elf %Plot phi component of two-way pattern 

meshc(Ur,Wr,Two_way_DAFp_db),grid,axis([-1 1-11 pmin,0]),grid 

axis square 

xlabel('U = sin\theta*cos\phi' ) 
ylabel('W = cosXtheta') 
zlabel (' Normalized F_\phi| (dB) ') 
view (45, 45) 

title(' Two-way Pattern') 

end 

%end case Calculate 

%Case calculate two-way gain of XMTR and RCVR arrays 
case 'Gain' 
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%Get user input for integration intervals for theta and phi 
prompt = { 'Number of intervals for theta: . 

'Number of intervals for phi:'}; 
dlg_title = 'Set number of integration points'; 
num_lines = 1; 
def = { ' 4 ' , '4' }; 

answer = inputdlg(prompt,dlg_title,num_lines,def); 

% set the number of integration intervals (nt points per interval) 
ndivt=str2num(answer{1}); 
ndivp=str2num(answer{2}); 


%Get parameters 

h_eltype = get(findobj(gcf, 'Tag ',' eltype '),' Value ') ; 


if 

h_eltype 

== 1, 

diptype = 

' h' ; 

: end 

if 

h_eltype 

== 2, 

diptype = 

's'; 

: end 

if 

h_eltype 

== 3, 

diptype = 

' i' J 

: end 


h_dipdir = get(findobj(gcf, 'Tag ',' dipdir '),' Value' ); 
if h_dipdir == 1, dipdir = 'x'; end 
if h_dipdir == 2, dipdir = 'z' ; end 
f = (str2num(get(findobj(gcf,'Tag','frequency'),'String')))*le6; 
h = (str2num(get(findobj(gcf,' Tag' , 'height' ), 'String '))) ; 
thetas = (str2num(get(findobj(gcf, 'Tag ',' thetas '),' String '))) ; 
phis = (str2num(get(findobj(gcf ,' Tag ',' phis '),' String '))) ; 


Nxt = str2num(get 
Nzt = str2num(get 
dxt = str2num(get 
dzt = str2num(get 
Mxt = str2num(get 
Mzt = str2num(get 
lxt = str2num(get 
lzt = str2num(get 

Nxr = str2num(get 
Nzr = str2num(get 
dxr = str2num(get 
dzr = str2num(get 
Mxr = str2num(get 
Mzr = str2num(get 
lxr = str2num(get 
lzr = str2num(get 


(findobj(gcf, 'Tag' , 
(findobj(gcf, 'Tag' , 
(findobj(gcf, 'Tag' , 
(findobj(gcf, 'Tag' , 
(findobj(gcf, 'Tag' , 
(findobj(gcf, 'Tag' , 
(findobj(gcf, 'Tag' , 
(findobj(gcf, 'Tag' , 

(findobj(gcf, 'Tag' , 
(findobj(gcf, 'Tag' , 
(findobj(gcf, 'Tag' , 
(findobj(gcf, 'Tag' , 
(findobj(gcf, 'Tag' , 
(findobj(gcf, 'Tag' , 
(findobj(gcf, 'Tag' , 
(findobj(gcf, 'Tag' , 


Nxt' ), 

'String 1 

); 

'Nzt' ), 

'String' 

); 

dxt' ) , 

'String' 

); 

dzt' ) , 

'String' 

); 

Mxt' ), 

'String' 

); 

Mzt ' ) , 

'String' 

); 

'lxt' ), 

'String' 

); 

lzt' ), 

'String' 

); 

’Nxr ' ), 

'String' 

); 

’Nzr ' ), 

'String' 

); 

dxr' ), 

'String' 

); 

dzr ' ), 

'String' 

); 

’Mxr' ), 

'String' 

); 

'Mzr' ), 

'String' 

); 

lxr'), 

'String' 

); 

lzr' ), 

'String' 

); 


txs_xdist = get(findobj(gcf ,' Tag ',' txs_xdist '),' Value ’) ; 
txs_peddbxval = get(findobj(gcf, 'Tag ',' txs_peddbx' ), 'Value' ); 
if txs_peddbxval==l, txs_peddbx=15; end 

txs_peddbx=20; end 
txs_peddbx=25; end 
txs_peddbx=30; end 
txs_peddbx=35; end 
txs_peddbx=40; end 
txs_peddbx=45; end 
txs_peddbx=50; end 


if 

if 

if 

if 

if 

if 

if 


txs_peddbxval==2, 
txs_peddbxval==3, 
txs_peddbxval==4, 
txs_peddbxval==5, 
txs_peddbxval==6, 
txs_peddbxval==7, 
txs_peddbxval==8, 


txs_nexpx = get(findobj(gcf ,' Tag ',' txs_nexpx '),' Value ') ; 


txs_zdist = get(findobj(gcf, 'Tag ',' txs_zdist '),' Value ') ; 
txs_peddbzval = get(findobj(gcf, 'Tag ',' txs_peddbz '),' Value '); 
if txs_peddbzval==l, txs_peddbz=15; end 
if txs_peddbzval==2, txs_peddbz=20; end 
if txs_peddbzval==3, txs_peddbz=25; end 
if txs_peddbzval==4, txs_peddbz=30; end 
if txs_peddbzval==5, txs_peddbz=35; end 
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if txs_peddbzval==6, txs_peddbz=40; end 
if txs_peddbzval==7, txs_peddbz=45; end 
if txs_peddbzval==8, txs_peddbz=50; end 
txs_nexpz = get(findobj(gcf, 'Tagtxs_nexpzValue ') ; 


tx_xdist = get(findobj(gcf Tagtx_xdist '), 'Value ') ; 
tx_peddbxval = get(findobj(gcf ,' Tagtx_peddbx '),' Value ') ; 
if tx_peddbxval==l, tx_peddbx=15; end 
tx peddbx=20; end 
tx_peddbx=25; end 
tx_peddbx=30; end 
tx_peddbx=35; end 
tx_peddbx=40; end 
tx_peddbx=45; end 
tx_peddbx=50; end 


if tx_peddbxval==2, 
if tx_peddbxval==3, 
if tx_peddbxval==4, 
if tx_peddbxval==5, 
if tx_peddbxval==6, 
if tx_peddbxval==7, 
if tx_peddbxval==8, 


tx_nexpx = get(findobj(gcf Tagtx_nexpxValue ') ; 


tx_zdist = get(findobj(gcf Tagtx_zdist '), 'Value ') ; 
tx_peddbzval = get(findobj(gcf ,' Tagtx_peddbz '),' Value ') ; 
if tx_peddbzval==l, tx peddbz=15; end 
tx_peddbz=20; end 
tx_peddbz=25; end 
tx_peddbz=30; end 
tx_peddbz=35; end 
tx_peddbz=40; end 
tx_peddbz=45; end 
tx_peddbz=50; end 


if tx_peddbzval==2, 
if tx_peddbzval==3, 
if tx_peddbzval==4, 
if tx_peddbzval==5, 
if tx_peddbzval==6, 
if tx_peddbzval==7, 
if tx_peddbzval==8, 


tx_nexpz = get(findobj(gcf Tagtx_nexpz '), 'Value ') ; 


txs_ampx = qetamplitudes(Mxt,txs xdist.txs peddbx,txs nexpx) 
txs_ampz = getamplitudes(Mzt,txs_zdist,txs_peddbz,txs_nexpz) 
tx_ampx = getamplitudes(Nxt,tx_xdist,tx_peddbx,tx_nexpx); 
tx_ampz = getamplitudes(Nzt,tx_zdist,tx_peddbz,tx_nexpz); 


rxs_xdist = get(findobj(gcf Tagrxs_xdistValue ’) ; 
rxs_peddbxval = get(findobj(gcf, 'Tag'rxs_peddbx '), 'Value' ) 
if rxs_peddbxval==l, rxs_peddbx=15; end 

rxs_peddbx=20; end 
rxs_peddbx=25; end 
rxs_peddbx=30; end 
rxs_peddbx=35; end 
rxs_peddbx=40; end 
rxs_peddbx=45; end 
rxs_peddbx=50; end 


if 

if 

if 

if 

if 

if 

if 


rxs_peddbxval==2, 
rxs_peddbxval==3, 
rxs_peddbxval==4, 
rxs peddbxval==5, 
rxs peddbxval==6, 
rxs_peddbxval==7, 
rxs_peddbxval==8, 


rxs_nexpx = get(findobj(gcf ,' Tag ',' rxs_nexpx '),' Value ') ; 


rxs_zdist = get(findobj(gcf ,' Tag ',' rxs_zdist '),' Value ') ; 
rxs_peddbzval = get(findobj(gcf,' Tag ',' rxs_peddbz '),' Value' ) 


if 

rxs_peddbzval= 

= 1, 

rxs_peddbz=15; 

end 

if 

rxs_peddbzval= 

=2, 

rxs peddbz=20; 

end 

if 

rxs peddbzval= 

=3, 

rxs peddbz=25; 

end 

if 

rxs_peddbzval= 

= 4, 

rxs_peddbz=30; 

end 

if 

rxs_peddbzval= 

=5, 

rxs_peddbz=35 ; 

end 

if 

rxs_peddbzval= 

= 6, 

rxs_peddbz=40; 

end 

if 

rxs_peddbzval= 

=7, 

rxs_peddbz=45; 

end 

if 

rxs_peddbzval= 

= 8, 

rxs_peddbz=50; 

end 


rxs_nexpz = get(findobj(gcf ,' Tag ',' rxs_nexpz '),' Value ') ; 


rx_xdist = get(findobj(gcf ,' Tagrx_xdist '), 'Value '); 
rx_peddbxval = get(findobj(gcf ,' Tag ',' rx_peddbx '),' Value ') ; 
if rx_peddbxval==l, rx_peddbx=15; end 
if rx_peddbxval==2, rx_peddbx=20; end 
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if rx_peddbxval==3, rx_peddbx=25; end 
if rx_peddbxval==4, rx_peddbx=30; end 
if rx_peddbxval==5, rx_peddbx=35; end 
if rx_peddbxval==6, rx peddbx=40; end 
if rx_peddbxval==7, rx_peddbx=45; end 
if rx_peddbxval==8, rx_peddbx=50; end 
rx_nexpx = get(findobj(gcf, 'Tagrx_nexpx '), 'Value ') ; 

rx_zdist = get(findobj(gcf Tagrx_zdist '), 'Value ') ; 


rx_peddbzval = get(findobj(gcf ,' Tag' , 

'rx_peddbz 

if 

rx_peddbzval= 

= 1, 

rx_peddbz=15; 

end 

if 

rx_peddbzval= 

=2, 

rx_peddbz=20; 

end 

if 

rx_peddbzval= 

=3, 

rx_peddbz=25; 

end 

if 

rx_peddbzval= 

= 4, 

rx_peddbz=30 ; 

end 

if 

rx_peddbzval= 

= 5, 

rx_peddbz=35 ; 

end 

if 

rx_peddbzval= 

= 6, 

rx_peddbz=40; 

end 

if 

rx_peddbzval= 

= 7, 

rx_peddbz=45; 

end 

if 

rx_peddbzval= 

= 8, 

rx_peddbz=50; 

end 


rx_nexpz = get(findobj(gcf, 'Tag ',' rx_nexpz '), 'Value ') ; 

rxs_ampx = getamplitudes(Mxr,rxs_xdist,rxs_peddbx,rxs_nexpx); 
rxs_ampz = getamplitudes(Mzr,rxs zdist.rxs peddbz,rxs nexpz); 
rx_ampx = getamplitudes(Nxr,rx_xdist,rx_peddbx,rx_nexpx); 
rx_ampz = getamplitudes(Nzr,rx_zdist,rx_peddbz,rx_nexpz); 

% Compute Gain of Transmit DSA 

[txgain,txemax,txprad]=compute_gain(' t' ,f,dxt,Nxt, dzt,Nzt,lxt,Mxt,lzt,Mzt,h,dip 
type,dipdir,thetas,phis,txs_ampx,txs_ampz,tx_ampx,tx_ampz,ndivt,ndivp); 

% Compute Gain of Receive DSA 

[rxgain,rxemax,rxprad]=compute_gain('r',f,dxr,Nxr,dzr,Nzr,lxr,Mxr,lzr,Mzr,h,dip 
type,dipdir,thetas,phis,rxs_ampx,rxs_ampz,rx_ampx,rx_ampz,ndivt,ndivp); 

txgdb=10*logl0(txgain); 

disp ( [ 'total radiated power, prad = ',num2str(txprad)]) 

disp(['max field value of transmit array, V/m = ', num2str(txemax)]) 

disp ([' transmit gain = ',num2str(txgain),' , in dB = ' ,num2str(txgdb)]) 

rxgdb=10*logl0(rxgain); 

disp (['max field value of receive array, V/m = ', num2str(rxemax)]) 
disp ([' receive gain = ', num2str(txgain) ,' , in dB = ' ,num2str(rxgdb) ]) 

% Compute two-way gain 
tway_gain=txgain*rxgain; 
tway_gaindb=10*logl0(tway_gain); 

nl = sprintf ( '\n' ); 

msgl = 'Transmit Array : '; 

msg2 = ['Numeric gain = ',num2str(txgain) ,.. . 

', Gain in dB = ',num2str(txgdb),' dB' ]; 
msg3 = 'Receive Array 

msg4 = ['Numeric gain = ',num2str(rxgain) ,.. . 

', Gain in dB = ',num2str(rxgdb),' dB' ]; 
msg5 = ['Two-way numeric gain = ' ,num2str(tway_gain),nl,... 

'Two-way Gain in dB = ' ,num2str(tway_gaindb),' dB'] ; 

h = msgbox([msgl nl msg2 nl nl msg3 nl msg4 nl nl msg5],... 

'Gain Results' ); 

%end case Gain 
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case 'Close' 

h_figs = get(0, 'children ') ; 
for fig = h_figs' 
delete(fig); 

end 

end %switch 


^Functions for input validation^ 


% validates Frequency 

function freq_out = getFreq(freq) 


temp = str2num(freq); 
if (isempty(temp)) | (temp <= 0) 

errordlg ('Enter a postive number.', ... 

'Frequency setting', 'error'); 
temp = 100; % default set to 100 

end 

freq_out = temp; 

% end getFreq 

% validates phi starting angle 
function o_pstart = getPStart(start) 


tempi = str2num(start); 

temp2 = str2num(get(findobj(gcf, 'Tag' , 'pstop' ), 'String' )); 

if (isempty(tempi)) | tempi < 0 I tempi > 180 

errordlg ('Enter a Phi Starting angle between 0 and 180 degrees 
'Angle Status', 'error'); 
tempi =90; % default phi start angle 
elseif (start == 'i' | start == 'j') 

errordlg ('Enter a Phi Starting angle between 0 and 180 degrees 
'Angle Status', 'error'); 
tempi =0; % default phi start angle 
elseif tempi > temp2 % phi start greater than phi stop angle 

errordlg ('Phi starting angle is greater than ending angle!', 
'Angle Status ',' error ') ; 
tempi =0; % default phi starting angle 

end 

o_pstart = tempi; 

% end getPStart 

% validates phi ending angle 
function o_pstop = getPStop(stop) 


tempi = str2num(stop); 

temp2 = str2num(get(findobj(gcf ,' Tag ',' pstart '),' String ')); 


if (isempty(tempi)) | tempi < 0 I tempi > 180 

errordlg ('Enter a Phi ending angle between starting angle and 180 
degrees.' , ... 


'Angle Status', 'error'); 

tempi = 180; 

elseif temp2 > tempi % phi start greater than phi stop angle 
errordlg ('Phi ending angle is less than starting angle!', ... 

'Angle Status ',' error ') ; 
tempi = 180; % default phi ending angle 
elseif (stop == 'i' I stop == 'j') 

errordlg ('Enter a Phi ending angle between 0 and 180 degrees.', 
'Angle Status', 'error'); 
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tempi = 180; % default phi ending angle 

end 

o_pstop = tempi; 

% end getPStop 

% validates theta starting angle 
function o_tstart = getTStart(start) 


tempi = str2num(start); 

temp2 = str2num(get(findobj(gcf, 'TagtstopString ')) ; 


if (isempty(tempi)) | tempi < 0 I tempi > 180 

errordlg ('Enter a Theta Starting angle between 
'Angle Status', 'error'); 
tempi =0; % default theta start angle 
elseif (start == 'i' | start == 'j') 

errordlg ('Enter a Theta Starting angle between 
'Angle Status', 'error'); 
tempi = 0; % default theta start angle 
elseif tempi > temp2 % theta start greater than 
errordlg ('Theta starting angle is greater than 
'Angle Status','error'); 
tempi =0; % default theta ending angle 

end 

o_tstart = tempi; 

% end getTStart 


0 and 180 degrees.'. 


0 and 180 degrees.'. 


theta stop angle 
ending angle! ' , 


% validates theta ending angle 
function o_tstop = getTStop(stop) 


tempi = str2num(stop); 

temp2 = str2num(get(findobj(gcf ,' Tag ',' tstart '),' String ')) ; 

if (isempty(tempi)) | (tempi < 0) I (tempi > 180) 

errordlg (' Enter a Theta ending angle between starting angle and 180 
degrees.' , ... 

'Theta Stop Angle Status', 'error'); 
tempi = 180; % default theta stop angle 

elseif (temp2 > tempi) % theta stop less than theta start angle 
errordlg ('Theta ending angle is less than starting angle!', ... 

'Theta Stop Angle Status ',' error ') ; 
tempi = 180; % default theta stop angle 
elseif (stop == 'i' I stop == 'j') 

errordlg ('Enter a Theta Ending angle between 0 and 180 degrees.', 
'Theta Stop Angle Status', 'error'); 
tempi = 180; % default theta ending angle 

end 

o_tstop = tempi; 

% end getTStop 

% validates Number of Elements in array or subarray 
function nel_out = getNEL(el_str) 


temp = str2double(el_str); 
if (isempty(temp)) || (floor(temp) <= 0) 

errordlg ('Enter a postive integer.', ... 

'Number of Elements', 'error'); 
temp =2; % default set to 2 

end 

nel_out = floor(temp); 

% end getNEL 
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% validates Array Element Spacing 
function o_del = getSpacing(str) 

temp = str2double(str); 

if (isempty(temp)) || (floor(temp) < 0) 

errordlg ('Please check spacing of array elements.',.. 

' Element Spacing', 'error'); 
temp =1; % default 

end 

o_del = temp; 

% end getSpacing 

% validates Scan Angle for Theta and Phi 
function o_scan = getScanAngle(scan_angle) 

temp = scan_angle; 

if (isempty(temp)) | (temp < 0) I (temp > 180) 

errordlg (' Enter a Scan Angle between 0 and 180 degrees 
' Scan Angle', 'error'); 
temp = 90; % default scan angle 

end 

o_scan = temp; 

% end getScanAngle 
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function [SAFtheta,SAFphi,DAFtheta,DAFphi,U,W]=caf_hdip2(tr,f,dx, Nx, .. 
dz,Nz,lx,Mx,lz,Mz,h,dipdir,thetas,phis,tstart,tstop, ... 
pstart,pstop,delt,delp,ampxs,ampzs,ampx,ampz) 

% caf_hdip2.m 
% Version: 2.2 
% Author: Cher Hock Hin 
% Advisor: Professor David C. Jenn 
% Date: 16 August 2012 

% Function computes array pattern for DSA in phi and theta cut (non-dB 
% Linear array of half-way dipoles with a ground plane. 

% z is array axis; ground plane is the xz plane. 

% y is normal to ground plane. 

% Function inputs: 

% f = frequency 

% dx = Subarray element spacing in x-direction in wavelengths 
% Nx = Number of subarray elements in x-direction 
% dz = Subarray element spacing in z-direction in wavelengths 
% Nz = Number of subarray elements in z-direction 
% lx = Subarray spacing in x-direction in wavelengths 

% Mx = Number of subarrays in x-direction 

% lz = Subarray spacing in z-direction in wavelengths 

% Mz = Number of subarrays in z-direction 

% h = height of element above ground plane in wavelengths 
% dipdir = dipole direction ("x"=x-direction,"z"=z-direction) 

% thetas = theta scan angle 
% phis = phi scan angle 
% tstart = start angle for theta 
% tstop = stop angle for theta 
% pstart = start angle for phi 
% pstop = stop angle for phi 

% delt = theta step size for pattern calculation 
% delp = phi step size for pattern calculation 
% ampxs = DSA amplitude distribution for x-plane 
% ampzs = DSA amplitude distribution for z-plane 
% ampx = Subarray amplitude distribution for x-plane 
% ampz = Subarray amplitude distribution for z-plane 

% Function outputs: 

% SAFtheta = Array pattern of single subarray for theta (complex) 

% SAFphi = Array pattern of single subarray for phi (complex) 

% DAFtheta = Array pattern of DSA for theta (complex) 

% DAFphi = Array pattern of DSA for theta (complex) 

warning off 
rad=pi/180; 

wave=3e8/f; %find wavelength 
beta=2*pi/wave; %find beta 
bk = 2*pi; 

%Set for loop values 

it=floor ( (tstop-tstart)/delt)+1; 

ip=floor((pstop-pstart)/delp)+1; 

%0pen waitbar 

if tr == 't', txrx = 'transmit'; end 
if tr == 'r', txrx = 'receive'; end 

msg = ['Computing ' ,txrx, ' array pattern for ' ,num2str(it*ip), ' angles 
hwait=waitbar(0,msg); 
pause(0.1); 
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% determine the phase distribution 

us = sin(thetas*rad)*cos(phis*rad); 

ws = cos(thetas*rad); 

psix = bk*dx*us; 

psiz = bk*dz*ws; 

psix_s = bk*lx*us; 

psiz_s = bk*lz*ws; 

% generate exact phase required at each element. 

% positive scan corresponds to increasing phase lag with increasing n. 

xsix = -(2*[l:Nx] - (Nx + l))/2*psix; 

xminx = min(xsix (1),xsix (Nx)); 

xsix (1:Nx)=xsix(1:Nx)-xminx; 

qphx=xsix; 

xsiz = -(2*[l:Nz] - (Nz + l))/2*psiz; 
xminz = min(xsiz(1),xsiz(Nz)); 
xsiz(1:Nz)=xsiz(1:Nz)-xminz; 
qphz=xsiz; 

% generate exact phase required at each subarray. 

% positive scan corresponds to increasing phase lag with increasing n. 

xsix_s = -(2*[l:Mx] - (Mx + 1))/2*psix_s; 

xminx_s = min(xsix_s(1),xsix_s(Mx)); 

xsix_s(1:Mx)=xsix_s(1:Mx)-xminx_s; 

qphx_s=xsix_s; 

xsiz_s = —(2*[1:Mz] - (Mz + 1))/2*psiz_s; 
xminz_s = min(xsiz_s(1),xsiz_s(Mz)) ; 
xsiz_s (1:Mz) =xsiz_s (1:Mz) -xminz_s; 
qphz_s=xsiz_s; 


%Preallocate array size 

phi=zeros(ip,it); theta=zeros (ip,it); 

U=zeros (ip,it); V=zeros(ip,it); W=zeros(ip,it); 
AF=zeros (ip,it); AFs = zeros(ip,it); 
SAFtheta=zeros(ip,it); SAFphi=zeros(ip,it); 
DAFtheta=zeros(ip,it); DAFphi=zeros(ip,it); 
count=0; 

% begin Pattern loop 
for il=l:ip %Loop phi 

for i2=l:it %Loop theta 
figure (hwait); 
count=count+l; 

waitbar(count/(ip*it),hwait); 
phi(il,i2) = pstart + (il - l)*delp; 
phr = phi(il,i2)*rad; 

theta(il,i2) = tstart + (12 - l)*delt; 
thr = theta(il,i2)*rad; 
st = sin (thr); ct = cos (thr); 

cp = cos (phr); sp = sin (phr); 

u = st*cp; v = st*sp; 


U (i 1, i 2) = u ; 


V(il, 12) = v; 


w = ct ; 

W(il,i2)=w; 


sumx = 0; 

%Sum to get array factor of subarray- begin subarray loop - 

for n = 1:Nx 

nn = (2*n - (Nx + l))/2; 
argx = bk*dx*u*nn; 
for m = 1 :Nz 

mm = (2*m - (Nz+l))/2; 

argz = bk*dz*w*mm; 

phase = qphx (n) + qphz(m); 

sumx = sumx + ampx(n)*ampz(m)*exp(li*(phase+argx+argz)) 


end 
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end % end subarray loop - 

sumx_s = 0; 

%Sum to get array factor of distributed subarrays- begin dsa loop 

for n = l:Mx 

nn = (2*n - (Mx + l))/2; 
argx_s = bk*lx*u*nn; 
for m = 1 :Mz 

mm = (2*m - (Mz+l))/2; 

argz_s = bk*lz*w*mm; 

phase_s = qphx_s(n) + qphz_s(m); 

sumx_s = sumx_s + 

ampxs(n)*ampzs(m)*exp(li*(phase_s+argx_s+argz_s)); 

end 

end % end dsa loop - 

AF(il,i2) = sumx; %Array factor for each subarray 
AFs(il,i2) = sumx_s; %Array factor for dsa 
GF = 2*sin (bk*h*v); %Ground plane factor 

if dipdir=='z' %Colinear dipoles (z directed) 
EFtheta=cos(pi*w/2)./sqrt(l-w. A 2+le-5); 

EFphi=0; 

end 

if dipdir=='x' %Parallel case (x directed) 

EFtheta=cp.*w.*cos (pi*u/2) ./ (l-u. A 2 + le-5); 

EFphi=-sp.*cos(pi*u/2) ./ (l-u. A 2 + le-5); 

end 

%Compute Subarray Pattern 

SAFtheta(il,i2) = AF(il,i2)*GF*EFtheta; 

SAFphi(il,i2) = AF(il,i2)*GF*EFphi; 

%Compute DSA Pattern 

DAFtheta(il,i2) = AF(il,i2)*GF*EFtheta*AFs(il,i2); 

DAFphi(il,i2) = AF(il,i2)*GF*EFphi*AFs(il,i2); 


end 

end % end of pattern loop 

% close waitbar 
close(hwait); 

return 
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function [SAFtheta,SAFphi,DAFtheta,DAFphi,U,W]=caf_sdip2(tr,f,dx, Nx, .. 
dz,Nz,lx,Mx,lz,Mz,h,dipdir,thetas,phis,tstart,tstop, ... 
pstart,pstop,delt,delp,ampxs,ampzs,ampx,ampz) 

% caf_sdip2.m 
% Version: 2.2 
% Author: Cher Hock Hin 
% Advisor: Professor David C. Jenn 
% Date: 16 August 2012 

% Function computes array pattern for DSA in phi and theta cut (non-dB 
% Linear array of short dipoles with a ground plane. 

% z is array axis; ground plane is the xz plane. 

% y is normal to ground plane. 

% Function inputs: 

% f = frequency 

% dx = Subarray element spacing in x-direction in wavelengths 
% Nx = Number of subarray elements in x-direction 
% dz = Subarray element spacing in z-direction in wavelengths 
% Nz = Number of subarray elements in z-direction 
% lx = Subarray spacing in x-direction in wavelengths 

% Mx = Number of subarrays in x-direction 

% lz = Subarray spacing in z-direction in wavelengths 

% Mz = Number of subarrays in z-direction 

% h = height of element above ground plane in wavelengths 
% dipdir = dipole direction ("x"=x-direction,"z"=z-direction) 

% thetas = theta scan angle 
% phis = phi scan angle 
% tstart = start angle for theta 
% tstop = stop angle for theta 
% pstart = start angle for phi 
% pstop = stop angle for phi 

% delt = theta step size for pattern calculation 
% delp = phi step size for pattern calculation 
% ampxs = DSA amplitude distribution for x-plane 
% ampzs = DSA amplitude distribution for z-plane 
% ampx = Subarray amplitude distribution for x-plane 
% ampz = Subarray amplitude distribution for z-plane 

% Function outputs: 

% SAFtheta = Array pattern of single subarray for theta (complex) 

% SAFphi = Array pattern of single subarray for phi (complex) 

% DAFtheta = Array pattern of DSA for theta (complex) 

% DAFphi = Array pattern of DSA for theta (complex) 

warning off 
rad=pi/180; 

wave=3e8/f; %find wavelength 
beta=2*pi/wave; %find beta 
bk = 2*pi; 

%Set for loop values 

it=floor ( (tstop-tstart)/delt)+1; 

ip=floor((pstop-pstart)/delp)+1; 

%0pen waitbar 

if tr == 't', txrx = 'transmit'; end 
if tr == 'r', txrx = 'receive'; end 

msg = ['Computing ' ,txrx, ' array pattern for ' ,num2str(it*ip), ' angles 
hwait=waitbar(0,msg); 
pause(0.1); 
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% determine the phase distribution 

us = sin(thetas*rad)*cos(phis*rad); 

ws = cos(thetas*rad); 

psix = bk*dx*us; 

psiz = bk*dz*ws; 

psix_s = bk*lx*us; 

psiz_s = bk*lz*ws; 

% generate exact phase required at each element. 

% positive scan corresponds to increasing phase lag with increasing n. 

xsix = -(2*[l:Nx] - (Nx + l))/2*psix; 

xminx = min(xsix (1),xsix (Nx)); 

xsix (1:Nx)=xsix(1:Nx)-xminx; 

qphx=xsix; 

xsiz = -(2*[l:Nz] - (Nz + l))/2*psiz; 
xminz = min(xsiz(1),xsiz(Nz)); 
xsiz(1:Nz)=xsiz(1:Nz)-xminz; 
qphz=xsiz; 

% generate exact phase required at each subarray. 

% positive scan corresponds to increasing phase lag with increasing n. 

xsix_s = -(2*[l:Mx] - (Mx + 1))/2*psix_s; 

xminx_s = min(xsix_s(1),xsix_s(Mx)); 

xsix_s(1:Mx)=xsix_s(1:Mx)-xminx_s; 

qphx_s=xsix_s; 

xsiz_s = —(2*[1:Mz] - (Mz + 1))/2*psiz_s; 
xminz_s = min(xsiz_s(1),xsiz_s(Mz)) ; 
xsiz_s (1:Mz) =xsiz_s (1:Mz) -xminz_s; 
qphz_s=xsiz_s; 


%Preallocate array size 

phi=zeros(ip,it); theta=zeros (ip,it); 

U=zeros (ip,it); V=zeros(ip,it); W=zeros(ip,it); 
AF=zeros (ip,it); AFs = zeros(ip,it); 
SAFtheta=zeros(ip,it); SAFphi=zeros(ip,it); 
DAFtheta=zeros(ip,it); DAFphi=zeros(ip,it); 
count=0; 

% begin Pattern loop 
for il=l:ip %Loop phi 

for i2=l:it %Loop theta 
figure (hwait); 
count=count+l; 

waitbar(count/(ip*it),hwait); 
phi(il,i2) = pstart + (il - l)*delp; 
phr = phi(il,i2)*rad; 

theta(il,i2) = tstart + (12 - l)*delt; 
thr = theta(il,i2)*rad; 
st = sin (thr); ct = cos (thr); 

cp = cos (phr); sp = sin (phr); 

u = st*cp; v = st*sp; 


U (i 1, i 2) = u ; 


V(il, 12) = v; 


w = ct ; 

W(il,i2)=w; 


sumx = 0; 

%Sum to get array factor of subarray- begin subarray loop - 

for n = 1:Nx 

nn = (2*n - (Nx + l))/2; 
argx = bk*dx*u*nn; 
for m = 1 :Nz 

mm = (2*m - (Nz+l))/2; 

argz = bk*dz*w*mm; 

phase = qphx (n) + qphz(m); 

sumx = sumx + ampx(n)*ampz(m)*exp(li*(phase+argx+argz)) 


end 
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end % end subarray loop - 

sumx_s = 0; 

%Sum to get array factor of distributed subarrays- begin dsa loop 

for n = l:Mx 

nn = (2*n - (Mx + l))/2; 
argx_s = bk*lx*u*nn; 
for m = 1 :Mz 

mm = (2*m - (Mz+l))/2; 

argz_s = bk*lz*w*mm; 

phase_s = qphx_s(n) + qphz_s(m); 

sumx_s = sumx_s + 

ampxs(n)*ampzs(m)*exp(li*(phase_s+argx_s+argz_s)); 

end 

end % end dsa loop - 

AF(il,i2) = sumx; %Array factor for each subarray 
AFs(il,i2) = sumx_s; %Array factor for dsa 
GF = 2*sin (bk*h*v); %Ground plane factor 

if dipdir=='z' %Colinear dipoles (z directed) 

EFtheta=sqrt(1-w. A 2); 

EFphi=0; 

end 

if dipdir=='x' %Parallel case (x directed) 

EFtheta=cp.*w; 

EFphi=-sp; 

end 

%Compute Subarray Pattern 

SAFtheta(il,i2) = AF(il,i2)*GF*EFtheta; 

SAFphi(il,i2) = AF(il,i2)*GF*EFphi; 

%Compute DSA Pattern 

DAFtheta(il,i2) = AF(il,i2)*GF*EFtheta*AFs(il,i2); 

DAFphi(il,i2) = AF(il,i2)*GF*EFphi*AFs(il,i2); 


end 

end % end of pattern loop 

% close waitbar 
close(hwait); 

return 
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function [SAFtheta,SAFphi,DAFtheta,DAFphi,U,W]=caf_iso2(tr,f,dx,Nx, ... 
dz,Nz,lx,Mx,lz,Mz,h,thetas, phis,tstart,tstop, ... 
pstart,pstop,delt,delp,ampxs,ampzs,ampx,ampz) 

% caf_iso2.m 
% Version: 2.2 
% Author: Cher Hock Hin 
% Advisor: Professor David C. Jenn 
% Date: 16 August 2012 

% Function computes array pattern for DSA in phi and theta cut (non-dB). 
% Linear array of isotropic elements with a ground plane. 

% z is array axis; ground plane is the xz plane. 

% y is normal to ground plane. 

% Function inputs: 

% f = frequency 

% dx = Subarray element spacing in x-direction in wavelengths 
% Nx = Number of subarray elements in x-direction 
% dz = Subarray element spacing in z-direction in wavelengths 
% Nz = Number of subarray elements in z-direction 
% lx = Subarray spacing in x-direction in wavelengths 

% Mx = Number of subarrays in x-direction 

% lz = Subarray spacing in z-direction in wavelengths 

% Mz = Number of subarrays in z-direction 

% h = height of element above ground plane in wavelengths 
% thetas = theta scan angle 
% phis = phi scan angle 
% tstart = start angle for theta 
% tstop = stop angle for theta 
% pstart = start angle for phi 
% pstop = stop angle for phi 

% delt = theta step size for pattern calculation 
% delp = phi step size for pattern calculation 
% ampxs = DSA amplitude distribution for x-plane 
% ampzs = DSA amplitude distribution for z-plane 
% ampx = Subarray amplitude distribution for x-plane 
% ampz = Subarray amplitude distribution for z-plane 

% Function outputs: 

% SAFtheta = Array pattern of single subarray for theta (complex) 

% SAFphi = Array pattern of single subarray for phi (complex) 

% DAFtheta = Array pattern of DSA for theta (complex) 

% DAFphi = Array pattern of DSA for theta (complex) 

warning off 
rad=pi/180; 

wave=3e8/f; %find wavelength 
beta=2*pi/wave; %find beta 
bk = 2*pi; 

%Set for loop values 

it=floor((tstop-tstart)/delt)+1; 

ip=floor ( (pstop-pstart)/delp)+1; 

%0pen waitbar 

if tr == 't', txrx = 'transmit'; end 
if tr == 'r', txrx = 'receive'; end 

msg = ['Computing ' ,txrx, ' array pattern for ' ,num2str(it*ip), ' angles'] 
hwait=waitbar(0,msg); 
pause(0.1); 

% determine the phase distribution 
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us = sin(thetas*rad)*cos(phis*rad); 

ws = cos (thetas*rad); 

psix = bk*dx*us; 

psiz = bk*dz*ws; 

psix_s = bk*lx*us; 

psiz_s = bk*lz*ws; 

% generate exact phase required at each element. 

% positive scan corresponds to increasing phase lag with increasing n. 

xsix = -(2*[l:Nx] - (Nx + l))/2*psix; 

xrainx = min(xsix(1),xsix (Nx)); 

xsix(1:Nx)=xsix(1:Nx)-xminx; 

qphx=xsix; 

xsiz = ~(2*[l:Nz] - (Nz + l))/2*psiz; 
xminz = min(xsiz(1),xsiz(Nz)); 
xsiz(1:Nz)=xsiz(1:Nz)-xminz; 
qphz=xsiz; 

% generate exact phase required at each subarray. 

% positive scan corresponds to increasing phase lag with increasing n. 

xsix_s = ~(2*[l:Mx] - (Mx + 1))/2*psix_s; 

xminx_s = min(xsix_s(1),xsix_s(Mx)); 

xsix_s(1:Mx)=xsix_s(1:Mx)-xminx_s; 

qphx_s=xsix_s; 

xsiz_s = —(2 *[1:Mz] - (Mz + 1))/2*psiz_s; 
xminz_s = min(xsiz_s(1),xsiz_s(Mz)); 
xsiz_s(1:Mz)=xsiz_s(1:Mz)-xminz_s; 
qphz_s=xsiz_s; 


%Preallocate array size 

phi=zeros(ip,it); theta=zeros (ip,it); 

U=zeros (ip,it); V=zeros(ip,it); W=zeros(ip,it); 
AF=zeros (ip,it); AFs = zeros(ip,it); 
SAFtheta=zeros(ip,it); SAFphi=zeros(ip,it); 
DAFtheta=zeros(ip,it); DAFphi=zeros(ip,it); 
count=0; 

% begin Pattern loop 
for il=l:ip %Loop phi 

for i2=l:it %Loop theta 
figure (hwait); 
count=count+l; 

waitbar(count/(ip*it),hwait); 
phi(il,i2) = pstart + (il - l)*delp; 
phr = phi(il,i2)*rad; 

theta(il,i2) = tstart + (i2 - l)*delt; 
thr = theta(il,i2)*rad; 
st = sin (thr); ct = cos (thr); 

cp = cos (phr) ; sp = sin (phr) ; 

u = st*cp; v = st*sp; 


U (i 1, i 2) = u ; 


V(il, i2) = v; 


w = ct; 

W (i 1, i 2) = w ; 


sumx = 0; 

%Sum to get array factor of subarray- begin subarray loop - 

for n = 1:Nx 

nn = (2*n - (Nx + l))/2; 
argx = bk*dx*u*nn; 
for m = 1 :Nz 

mm = (2*m - (Nz+l))/2; 

argz = bk*dz*w*mm; 

phase = qphx(n) + qphz(m); 

sumx = sumx + ampx(n)*ampz(m)*exp(li*(phase+argx+argz)) 

end 


end % end subarray loop 
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sumx_s = 0; 

%Sum to get array factor of distributed subarrays- begin dsa loop 

for n = l:Mx 

nn = (2*n - (Mx + l))/2; 
argx_s = bk*lx*u*nn; 
for m = 1 :Mz 

mm = (2*m - (Mz+l))/2; 

argz_s = bk*lz*w*mm; 

phase_s = qphx_s(n) + qphz_s(m); 

sumx_s = sumx_s + 

ampxs(n)*arapzs(m)*exp(1i*(phase_stargx_s+argz_s)); 

end 

end % end dsa loop - 

AF(il,i2) = sumx; %Array factor for each subarray 
AFs(il,i2) = sumx_s; %Array factor for dsa 
GF = 2*sin (bk*h*v); %Ground plane factor 

%Compute Subarray Pattern 
SAFtheta(il,i2) = AF(il,i2)*GF; 

SAFphi(i1,i2) = AF (il, i2) *GF; 

%Compute DSA Pattern 

DAFtheta(il,i2) = AF(il,i2)*GF*AFs(il,i2); 

DAFphi(il,i2) = AF(il,i2)*GF*AFs(il,i2); 


end 

end % end of pattern loop 

% close waitbar 
close(hwait); 

return 
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function [arapnorm] = getamplitudes(N,dist,pedval,nexp) 

% filename: getamplitudes . m 

% Description: This program calculates the amplitude distribution. 

% Author: Prof. David C. Jenn 
% Modified by: Cher Hock Hin 
% Date: 4 August 2012 

% Required subroutines: tayl.m; cosine. m; bayliss.m 

% uniform array excitation coefficients (=l/nel) 
if dist == 1 

for i = 1:N 
amp(i) = 1/N; 
end 

end 

% call subroutine to compute taylor coefficients 
% NEL MUST BE EVEN 
if dist == 2 

amp(l:N) = tayl (N,pedval,nexp); 

end 

% call subroutine to compute cosine-on-a-pedestal distribution 
if dist == 3 

amp(l:N) = cosine(N,pedval,nexp); 

end 

% call subroutine to compute bayliss distribution for difference beams 
% NEL MUST BE EVEN 
if dist == 4 

amp(l:N) = bayliss(N,pedval,nexp); 

end 

% compute triangular coefficients 
if dist == 5 

for i - i :N 

n = (2*i - (N + 1) ) /2; 
amp(i) = 1 - abs (2*n/N); 

end 

end 

% normalize all coefficients to the maximum value 
ampmax = max(amp); 
ampnorm = amp/ampmax; 
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function amp=bayliss(nel,sll, nbar) 

% SUBROUTINE TO COMPUTE LINEAR BAYLISS COEFFICIENTS FOR DIFFERENCE 
% BEAMS. LINEAR BAYLISS DISTRIBUTION 

nbarl=nbar-l; 
for i=l:11 

mu(i)=(i-1)+.5; 
end 

sll25=sll-25; 
z (l)=1.87 + sll25*.038; 
z(2)=2.50+S1125*.016; 
z(3)=3.35+sll25*.019; 
z(4)=4.25+sll25*.016; 
a=l.45+S1125*.042; 
for ns=5:nbar 

z (ns)=sqrt(a A 2+ns A 2); 
end 

sigma=mu(nbar+1)/z(nbar); 
del=2/(nel-1); 
for i=l:nel/2 
rho=del/2+del*(i-1); 
for mms=l:nbar 
bb=l; 

for ns=l:nbarl 

bb=bb*(1-(mu(mms)/(sigma*z(ns))) A 2); 
end 
bbb=l; 

for lls=l:nbar 
if 11s ~= mms 

bbb=bbb*(1-(mu(mms)/mu(11s)) A 2); 
end 
end 

bes= (-1) A mms; 

b(mms)=mu(mms) A 2/bes*bb/bbb; 

end 

gg=0; 

for lls=l:nbar 

pmu=mu(11s)*pi*rho; 
bes=sin(pmu); 
gg=gg+b(11s)*bes; 
end 

gg=abs (gg); 
amp(nel/2+i)=gg; 
amp(nel/2+l-i)=gg; 
end 


function amp=cosine(nel, peddb, nexp) 

% cosine on a pedestal distribution 

% scale the number of elements to correspond to -1 through +1 

ped=10 A (-peddb/20.); 
for n=l:nel 

xn= (n—1)/ (nel-1)*2-1; 

amp(n)=(1.-ped)*abs(cos(xn*pi/2)) A nexp+ped; 

end 
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function amp=tayl(noel,slldb, nbar) 

% COMPUTES TAYLOR DISTRIBUTION FOR GIVEN SIDELOBE LEVEL AND NBAR 

% AMP=ARRAY OF AMPLITUDES COMPUTED BY SUBROUTINE 

% NOEL=NUMBER OF ARRAY ELEMENTS 

% SLL=SIDELOBE LEVEL IN DB 

% NBAR=NBAR IN TAYLOR DISTRIBUTION (.LE.50) 

for i=l:noel 
amp(i)=.5; 
end 

if nbar~=l % if nbar > 1 
dbamp=2 0/log(10); 
sll=exp(abs(slldb)/dbamp); 
as=log (sll + sqrt(sll A 2-l))/pi; 
as=as A 2; 

s=nbar A 2/(as+nbar A 2-nbar+.25); 
nll=nbar-l; 
for ii=l:nll 
al=ii A 2/s; 

f(ii)=1; 

for jj=l:nll 

f(ii)=f(ii)*(1-al/(as + jj A 2-j j+.25)); 
end 

for jj=l:ii 

f (ii)=f(ii)/(1 + ii/(nbar-jj) ) ; 
end 
end 

m2=noel/2; % noel assumed positive 
if 2*m2 <= noel 
dum=.5; 
for ii=l:nll 
dum=dum+f(ii); 
end 

amp(m2+l)=2*dum; 
end 

for ii=l:m2 
k=noel+l-ii; 
for jj=l:nll 

amp (ii) =amp(ii) +f(jj)*cos(pi* j j *(k-ii)/noel); 
end 

amp (ii) =2*amp (ii) ; 
amp(k)=amp(ii); 
end 
end 
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function 

[gain,prad,emax]=compute_gain(tr,f,dx,Nx,dz,Nz,lx,Mx,Iz,Mz,h,diptype,dipdir, the 
tas,phis,ampxs,ampzs,ampx,ampz,ndivt,ndivp) 

% compute_gain.m 
% Version: 1.1 
% Author: Cher Hock Hin 
% Advisor: Professor David C. Jenn 
% Date: 11 August 2012 

% Function computes gain for DSA in x-z plane with y normal. 

% Uses 20 points Gaussian quadrature integration. 

% Reads integration constants from gausq20.m 

% Function inputs: 

% tr = •t' for XMTR and 'r' for RCVR 
% f = frequency 

% dx = Subarray element spacing in x-direction in wavelengths 
% Nx = Number of subarray elements in x-direction 
% dz = Subarray element spacing in z-direction in wavelengths 
% Nz = Number of subarray elements in z-direction 
% lx = Subarray spacing in x-direction in wavelengths 

% Mx = Number of subarrays in x-direction 

% lz = Subarray spacing in z-direction in wavelengths 

% Mz = Number of subarrays in z-direction 

% h = height of element above ground plane in wavelengths 
% diptype = element type ("s"=short,"h"=half-wave,"i"=isotropic) 

% dipdir = dipole direction ("x"=x-direction,"z"=z-direction) 

% thetas = theta scan angle 
% phis = phi scan angle 

% ndivt = number of integration intervals for theta 
% ndivp = number of integration intervals for phi 
% ampxs = DSA amplitude distribution for x-plane 
% ampzs = DSA amplitude distribution for z-plane 
% ampx = Subarray amplitude distribution for x-plane 
% ampz = Subarray amplitude distribution for z-plane 

% Function outputs: 

% gain = Numerical gain of DSA 
% prad = Total radiated power of DSA 
% emax = Maximum field value of DSA 

rad=pi/180; 

wave=3e8/f; %find wavelength 
beta=2*pi/wave; %find beta 
bk = 2*pi; 

% load data for 20-point Gaussian quadrature integration 

load gausq20.m 

xt=gausq20(:,1); 

at=gausq20(:,2); 

nt=length(xt); 

Timel=cputime; 

% integration interval in theta (degrees) 

Sl=0*rad; 

S2=180*rad; 

% integration interval in phi (degrees) 

Ql=0*rad; 

Q2=180*rad; 

% generate integration points in theta and phi 
% ndivt and ndivp points will be used but the contribution at 
% ndivt+1 and ndivp+1 
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ds=(S2-Sl)/ndivt; 
for i=l:ndivt+l 

SS(i)=(i—1)*ds; 

disp([ 'i,SS(i)= ' ,num2str(i), ' , ' ,num2str(SS (i)) ]) 

end 

dq=(Q2-Q1)/ndivp; 
for i=l:ndivp+l 

QQ (i) = (i — 1)*dq; 

disp([ 'i, QQ(i) = ' ,num2str (i), ' , ' ,num2str(QQ (i)) ]) 

end 

% subintervals in phi 
nphi = 0; 

for ii=l:ndivp 
Pl=dq/2; 

P2=(QQ(ii+1)+QQ(ii))/2; 
for n=l:nt 

nphi=nphi+l; 

wphi(nphi)=at (n) ; 

phi(nphi)=Pl*xt(n)+P2; 

end 

end 

% subintervals in theta 
ntheta=0; 
for ii=l:ndivt 
Tl=ds/2; 

T2=(SS(ii+1)+SS(ii))/2; 
for i=l:nt 

ntheta=ntheta+l; 
wtheta(ntheta)=at(i); 
theta (ntheta)=Tl*xt(i)+T2; 

end 

end 

ninteg=ntheta*nphi; 
if tr=='t', tr= 'transmit' ; end 
if tr== 'r 1 , tr= 'receive' ; end 

msg=[ 'Computing gain for ' ,tr, 1 array' ,sprintf ( ' \n' ), . . . 

'over ' ,num2str(ninteg), 1 integration points']; 
hwait=waitbar(0,msg); 

% determine the phase distribution 
us = sin(thetas*rad)*cos(phis*rad); 
ws = cos (thetas*rad); 
psix = bk*dx*us; 
psiz = bk*dz*ws; 

% generate exact phase required at each element. 

% positive scan corresponds to increasing phase lag with increasing n. 

xsix = -(2*[l:Nx] - (Nx + l))/2*psix; 

xminx = min(xsix (1),xsix (Nx)); 

xsix (1:Nx)=xsix(1:Nx)-xminx; 

qphx=xsix; 

xsiz = -(2*[l:Nz] - (Nz + l))/2*psiz; 
xminz = min(xsiz(1),xsiz(Nz)); 
xsiz(1:Nz)=xsiz(1:Nz)-xminz; 
qphz=xsiz; 

% generate exact phase required at each subarray. 

% positive scan corresponds to increasing phase lag with increasing n. 
xsix_s = -(2*[l:Mx] - (Mx + l))/2*psix; 
xminx_s = min(xsix_s(1),xsix_s(Mx)); 
xsix_s(1:Mx)=xsix_s(1:Mx)-xminx_s; 
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qphx_s=xsix_s; 

xsiz_s = ~(2*[l:Mz] - (Mz + l))/2*psiz; 
xminz_s = min(xsiz_s(1),xsiz_s(Mz)); 
xsiz_s(1:Mz)=xsiz_s(1:Mz)-xminz_s; 
qphz_s=xsiz_s; 

% compute field at the integration points 
emax=0; sum_pwr=0; ism=0; 

% begin Pattern loop 
for iphi=l:nphi %Loop phi 

for itheta=l:ntheta %Loop theta 
ism=ism+l; 

waitbar(ism/ninteg,hwait); 

thr = theta(itheta); phr = phi(iphi); 

st = sin (thr); ct = cos (thr); 

cp = cos (phr); sp = sin (phr); 

u = st*cp; v = st*sp; w = ct; 

sumx = 0; 

%Sum to get array factor of subarray- begin subarray loop - 

for n = 1:Nx 

nn = (2*n - (Nx + 1)) /2; 
argx = bk*dx*u*nn; 
for m = 1 :Nz 

mm = (2*m - (Nz+l))/2; 

argz = bk*dz*w*mm; 

phase = qphx(n) + qphz(m); 

sumx = sumx + ampx(n)*ampz(m)*exp(li*(phase+argx+argz)) 

end 

end % end subarray loop - 

sumx_s = 0; 

%Sum to get array factor of distributed subarrays- begin dsa loop - 

for n = 1:Mx 

nn = (2*n - (Mx + 1))/2; 
argx_s = bk*lx*u*nn; 
for m = 1 :Mz 

mm = (2*m - (Mz+l))/2; 

argz_s = bk*lz*w*mm; 

phase_s = qphx_s(n) + qphz_s(m); 

sumx_s = sumx_s + 

ampxs(n)*ampzs(m)*exp(li*(phase_s+argx_s+argz_s)); 

end 

end % end dsa loop - 

GF = 2*sin (bk*h*v); %Ground plane factor 

if diptype=='i 1 %Compute element factor for isotropic 
EFtheta=l; EFphi=-l; 

end 

if diptype=='h' %Compute element factor for half-wave dipole 
if dipdir=='z' %Colinear dipoles (z directed) 
EFtheta=cos(pi*w/2)./sqrt(1-w. A 2+le-5); 

EFphi=0; 

end 

if dipdir=='x' %Parallel case (x directed) 

EFtheta=cp.*w.*cos(pi*u/2) ./ (l-u. A 2 + le-5); 

EFphi=-sp.*cos(pi*u/2) ./ (l-u. A 2 + le-5); 

end 

end 

if diptype=='s' %Compute element factor for short dipole 
if dipdir=='z' %Colinear dipoles (z directed) 
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EFtheta=sqrt(1-w. A 2); 

EFphi=0; 

end 

if dipdir=='x' %Parallel case (x directed) 
EFtheta=cp.*w; 

EFphi=-sp; 

end 

end 

% compute emag 

emagt = sumx*sumx_s*GF*EFtheta; 
emagp = sumx*sumx_s*GF*EFphi; 
emagsq = abs(emagt)+ abs(emagp) A 2; 

if emagsq > emax, emax=emagsq; end % keep track of emax 
xint = emagsq*st; 

sum_pwr = sum_pwr + wphi (iphi)*wtheta(itheta)*xint; 

end 

end % end of pattern loop 

% close waitbar 
close(hwait); 

prad=Tl*Pl*sum_pwr; % compute radiated power 
gain=4*pi*emax/prad; % compute gain 

Time2=cputime-Timel; 

disp ([' runtime for integration: ' ,num2str(Time2)]) 
return 
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function varargout = two wav pattern(varargin) 

% TWO_WAY_PATTERN M-file for two wav pattern.fig 

% TWO_WAY_PATTERN, by itself, creates a new TWO_WAY_PATTERN or raises the 

existing 

% singleton* . 


% H = TWO_WAY_PATTERN returns the handle to a new TWO_WAY_PATTERN or the 

handle to 

% the existing singleton*. 

% TWO_WAY_PATTERN ('CALLBACK', hObject,eventData,handles, ...) calls the 

local 

% function named CALLBACK in TWO_WAY_PATTERN . M with the given input 

arguments . 

% TWO_WAY_PATTERN( 'Property','Value',...) creates a new TWO_WAY_PATTERN or 

raises the 

% existing singleton*. Starting from the left, property value pairs are 

% applied to the GUI before two way pattern OpeninqFcn gets called. An 

% unrecognized property name or invalid value makes property application 

% stop. All inputs are passed to two way pattern OpeninqFcn via varargin. 

% *See GUI Options on GUIDE'S Tools menu. Choose "GUI allows only one 

% instance to run (singleton) ". 

% See also: GUIDE, GUIDATA, GUIHANDLES 


% Edit the above text to modify the response to help two way pattern 
% Last Modified by GUIDE v2 . 5 04-Aug-2012 21:11:53 


% Begin initialization code - DO NOT EDIT 
gui_Singleton = 1; 
gui_State = struct (' gui_Name' , 

'gui_Singleton' , 

'gui_OpeningFcn' , 

'gui_OutputFcn' , 

'gui_LayoutFcn' , 

'gui_Callback' , 
if nargin && ischar(varargin{1}) 

gui_State.gui_Callback = str2func(varargin{1}); 

end 


mfilename, ... 
gui_Singleton, ... 
@two_way_pattern_OpeningFcn, 
@two_way_pattern_OutputFcn, 

[],... 

[]) ; 


if nargout 

[varargout{1:nargout}] = gui_mainfcn(gui_State, varargin{:}); 

else 

gui_mainfcn(gui_State, varargin):}); 

end 

% End initialization code - DO NOT EDIT 


% - Executes just before two wav pattern is made visible. 

function two wav pattern OpeninqFcn(hObject, eventdata, handles, varargin) 
% This function has no output args, see OutputFcn. 

% hObject handle to figure 

% eventdata reserved - to be defined in a future version of MATLAB 
% handles structure with handles and user data (see GUIDATA) 

% varargin command line arguments to two_way_pattern (see VARARGIN) 


% Choose default command line output for two wav pattern 
handles.output = hObject; 


% Update handles structure 
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guidata(hObject, handles); 


% UIWAIT makes two wav pattern wait for user response (see UIRESUME) 
% uiwait(handles .figurel); 


% - Outputs from this function are returned to the command line. 

function varargout = two wav pattern OutputFcn (hObject, eventdata, handles) 
% varargout cell array for returning output args (see VARARGOUT); 

% hObject handle to figure 

% eventdata reserved - to be defined in a future version of MATLAB 
% handles structure with handles and user data (see GUIDATA) 


% Get default command line output from handles structure 
varargout{1} = handles.output; 


% - Executes on button press in pushbuttonl . 

function pushbuttonl_Callback(hObject, eventdata, 
% hObject handle to pushbuttonl (see GCBO) 

% eventdata reserved - to be defined in a future 
% handles structure with handles and user data 


handles) 

version of MATLAB 
(see GUIDATA) 


% - Executes on button press in pushbutton2 . 

function pushbutton2_Callback(hObject, eventdata, 
% hObject handle to pushbutton2 (see GCBO) 

% eventdata reserved - to be defined in a future 
% handles structure with handles and user data 


handles) 

version of MATLAB 
(see GUIDATA) 


% - Executes on selection change in eltype . 

function eltype_Callback(hObject, eventdata, handles) 

% hObject handle to eltype (see GCBO) 

% eventdata reserved - to be defined in a future version of MATLAB 
% handles structure with handles and user data (see GUIDATA) 


% Hints: contents = cellstr (get(hObject,' String ')) returns eltype contents as 
cell array 

% contents { get(hObject , 1 Value ') } returns selected item from eltype 

% - Executes during object creation, after setting all properties. 

function eltype_CreateFcn(hObject, eventdata, handles) 

% hObject handle to eltype (see GCBO) 

% eventdata reserved - to be defined in a future version of MATLAB 
% handles empty - handles not created until after all CreateFcns called 

% Hint: popupmenu controls usually have a white background on Windows. 

% See ISPC and COMPUTER. 

if ispc && isequal(get(hObject,'BackgroundColor'), 
get (0, ' defaultUicontrolBackgroundColor' )) 

set(hObject,' BackgroundColor' , 'white' ); 

end 


function dzt_Callback(hObject, eventdata, handles) 

% hObject handle to dzt (see GCBO) 

% eventdata reserved - to be defined in a future version of MATLAB 
% handles structure with handles and user data (see GUIDATA) 

% Hints: get (hObject, 'String') returns contents of dzt as text 

% str2double(get (hObject,' String ')) returns contents of dzt as a double 

% - Executes during object creation, after setting all properties. 

function dzt_CreateFcn(hObject, eventdata, handles) 

% hObject handle to dzt (see GCBO) 

% eventdata reserved - to be defined in a future version of MATLAB 
% handles empty - handles not created until after all CreateFcns called 
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% Hint: edit controls usually have a white background on Windows. 

% See ISPC and COMPUTER. 

if ispc && isequal(get(hObject,'BackgroundColor'), 
get (0, ' defaultUicontrolBackgroundColor' )) 

set(hObject,' BackgroundColor' , 'white' ); 

end 

function dxt_Callback(hObject, eventdata, handles) 

% hObject handle to dxt (see GCBO) 

% eventdata reserved - to be defined in a future version of MATLAB 
% handles structure with handles and user data (see GUIDATA) 

% Hints: get (hObject, 'String') returns contents of dxt as text 

% str2double(get(hObject,'String')) returns contents of dxt as a double 

% - Executes during object creation, after setting all properties. 

function dxt_CreateFcn(hObject, eventdata, handles) 

% hObject handle to dxt (see GCBO) 

% eventdata reserved - to be defined in a future version of MATLAB 
% handles empty - handles not created until after all CreateFcns called 

% Hint: edit controls usually have a white background on Windows. 

% See ISPC and COMPUTER. 

if ispc && isequal(get(hObject,'BackgroundColor'), 
get (0, ' defaultUicontrolBackgroundColor' )) 

set(hObject,' BackgroundColor' , 'white' ); 

end 

function Nzt_Callback(hObject, eventdata, handles) 

% hObject handle to Nzt (see GCBO) 

% eventdata reserved - to be defined in a future version of MATLAB 
% handles structure with handles and user data (see GUIDATA) 

% Hints: get (hObject, 'String') returns contents of Nzt as text 

% str2double(get(hObject,'String')) returns contents of Nzt as a double 

% - Executes during object creation, after setting all properties. 

function Nzt_CreateFcn(hObject, eventdata, handles) 

% hObject handle to Nzt (see GCBO) 

% eventdata reserved - to be defined in a future version of MATLAB 
% handles empty - handles not created until after all CreateFcns called 

% Hint: edit controls usually have a white background on Windows. 

% See ISPC and COMPUTER. 

if ispc && isequal(get(hObject,'BackgroundColor'), 
get (0, ' defaultUicontrolBackgroundColor' )) 

set(hObject,' BackgroundColor' , 'white' ); 

end 

function Nxt_Callback(hObject, eventdata, handles) 

% hObject handle to Nxt (see GCBO) 

% eventdata reserved - to be defined in a future version of MATLAB 
% handles structure with handles and user data (see GUIDATA) 

% Hints: get (hObject, 'String') returns contents of Nxt as text 

% str2double(get(hObject,'String')) returns contents of Nxt as a double 

% - Executes during object creation, after setting all properties. 

function Nxt_CreateFcn(hObject, eventdata, handles) 

% hObject handle to Nxt (see GCBO) 

% eventdata reserved - to be defined in a future version of MATLAB 
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% handles empty - handles not created until after all CreateFcns called 

% Hint: edit controls usually have a white background on Windows. 

% See ISPC and COMPUTER. 

if ispc && isequal(get(hObject,'BackgroundColor'), 
get (0, ' defaultUicontrolBackgroundColor' )) 

set(hObject,' BackgroundColor' , 'white' ); 

end 

function Mxt_Callback(hObject, eventdata, handles) 

% hObject handle to Mxt (see GCBO) 

% eventdata reserved - to be defined in a future version of MATLAB 
% handles structure with handles and user data (see GUIDATA) 

% Hints: get(hObject,'String') returns contents of Mxt as text 

% str2double(get(hObject,'String')) returns contents of Mxt as a double 

% - Executes during object creation, after setting all properties. 

function Mxt_CreateFcn(hObject, eventdata, handles) 

% hObject handle to Mxt (see GCBO) 

% eventdata reserved - to be defined in a future version of MATLAB 
% handles empty - handles not created until after all CreateFcns called 

% Hint: edit controls usually have a white background on Windows. 

% See ISPC and COMPUTER. 

if ispc && isequal(get(hObject,'BackgroundColor'), 
get (0, ' defaultUicontrolBackgroundColor' )) 

set(hObject,' BackgroundColor' , 'white' ); 

end 

function Mzt_Callback(hObject, eventdata, handles) 

% hObject handle to Mzt (see GCBO) 

% eventdata reserved - to be defined in a future version of MATLAB 
% handles structure with handles and user data (see GUIDATA) 

% Hints: get (hObject, 'String') returns contents of Mzt as text 

% str2double(get(hObject,'String')) returns contents of Mzt as a double 

% - Executes during object creation, after setting all properties. 

function Mzt_CreateFcn(hObject, eventdata, handles) 

% hObject handle to Mzt (see GCBO) 

% eventdata reserved - to be defined in a future version of MATLAB 
% handles empty - handles not created until after all CreateFcns called 

% Hint: edit controls usually have a white background on Windows. 

% See ISPC and COMPUTER. 

if ispc && isequal(get(hObject,'BackgroundColor'), 
get (0, ' defaultUicontrolBackgroundColor' )) 

set(hObject,' BackgroundColor' , 'white' ); 

end 

function lxt_Callback(hObject, eventdata, handles) 

% hObject handle to lxt (see GCBO) 

% eventdata reserved - to be defined in a future version of MATLAB 
% handles structure with handles and user data (see GUIDATA) 

% Hints: get (hObject, 'String') returns contents of lxt as text 

% str2double(get(hObject,'String')) returns contents of lxt as a double 

% - Executes during object creation, after setting all properties. 

function lxt_CreateFcn(hObject, eventdata, handles) 

% hObject handle to lxt (see GCBO) 
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% eventdata reserved - to be defined in a future version of MATLAB 
% handles empty - handles not created until after all CreateFcns called 

% Hint: edit controls usually have a white background on Windows. 

% See ISPC and COMPUTER. 

if ispc && isequal(get(hObject, 'BackgroundColor' ), 
get ( 0, 'defaultUicontrolBackgroundColor' )) 

set(hObject,' BackgroundColor' , 'white' ); 

end 

function lzt_Callback(hObject, eventdata, handles) 

% hObject handle to lzt (see GCBO) 

% eventdata reserved - to be defined in a future version of MATLAB 
% handles structure with handles and user data (see GUIDATA) 

% Hints: get(hObject,'String') returns contents of lzt as text 

% str2double(get(hObject,'String')) returns contents of lzt as a double 

% - Executes during object creation, after setting all properties. 

function lzt_CreateFcn(hObject, eventdata, handles) 

% hObject handle to lzt (see GCBO) 

% eventdata reserved - to be defined in a future version of MATLAB 
% handles empty - handles not created until after all CreateFcns called 

% Hint: edit controls usually have a white background on Windows. 

% See ISPC and COMPUTER. 

if ispc && isequal(get(hObject ,' BackgroundColor ') , 
get (0, 'defaultUicontrolBackgroundColor' )) 

set (hObject, ' BackgroundColor' , 'white' ); 

end 

function Mxr_Callback(hObject, eventdata, handles) 

% hObject handle to Mxr (see GCBO) 

% eventdata reserved - to be defined in a future version of MATLAB 
% handles structure with handles and user data (see GUIDATA) 

% Hints: get (hObject, 'String') returns contents of Mxr as text 

% str2double(get(hObject,'String')) returns contents of Mxr as a double 

% - Executes during object creation, after setting all properties. 

function Mxr_CreateFcn(hObject, eventdata, handles) 

% hObject handle to Mxr (see GCBO) 

% eventdata reserved - to be defined in a future version of MATLAB 
% handles empty - handles not created until after all CreateFcns called 

% Hint: edit controls usually have a white background on Windows. 

% See ISPC and COMPUTER. 

if ispc && isequal(get(hObject ,' BackgroundColor ') , 
get ( 0, 'defaultUicontrolBackgroundColor' )) 

set (hObject, ' BackgroundColor' , 'white' ); 

end 

function Mzr_Callback(hObject, eventdata, handles) 

% hObject handle to Mzr (see GCBO) 

% eventdata reserved - to be defined in a future version of MATLAB 
% handles structure with handles and user data (see GUIDATA) 

% Hints: get (hObject, 'String') returns contents of Mzr as text 

% str2double(get(hObject,'String')) returns contents of Mzr as a double 

% - Executes during object creation, after setting all properties. 

function Mzr_CreateFcn(hObject, eventdata, handles) 
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% hObject handle to Mzr (see GCBO) 

% eventdata reserved - to be defined in a future version of MATLAB 
% handles empty - handles not created until after all CreateFcns called 

% Hint: edit controls usually have a white background on Windows. 

% See ISPC and COMPUTER. 

if ispc && isequal(get(hObject BackgroundColor ') , 
get (0, ' defaultUicontrolBackgroundColor' )) 

set (hObject, ' BackgroundColor' , 'white' ); 

end 

function lxr_Callback(hObject, eventdata, handles) 

% hObject handle to lxr (see GCBO) 

% eventdata reserved - to be defined in a future version of MATLAB 
% handles structure with handles and user data (see GUIDATA) 

% Hints: get (hObject, 'String') returns contents of lxr as text 

% str2double(get(hObject,'String')) returns contents of lxr as a double 

% - Executes during object creation, after setting all properties. 

function lxr_CreateFcn(hObject, eventdata, handles) 

% hObject handle to lxr (see GCBO) 

% eventdata reserved - to be defined in a future version of MATLAB 
% handles empty - handles not created until after all CreateFcns called 

% Hint: edit controls usually have a white background on Windows. 

% See ISPC and COMPUTER. 

if ispc && isequal(get(hObject ,' BackgroundColor ') , 
get (0, ' defaultUicontrolBackgroundColor' )) 

set(hObject,' BackgroundColor' , 'white' ); 

end 

function lzr_Callback(hObject, eventdata, handles) 

% hObject handle to lzr (see GCBO) 

% eventdata reserved - to be defined in a future version of MATLAB 
% handles structure with handles and user data (see GUIDATA) 

% Hints: get (hObject, 'String') returns contents of lzr as text 

% str2double(get(hObject,'String')) returns contents of lzr as a double 

% - Executes during object creation, after setting all properties. 

function lzr_CreateFcn(hObject, eventdata, handles) 

% hObject handle to lzr (see GCBO) 

% eventdata reserved - to be defined in a future version of MATLAB 
% handles empty - handles not created until after all CreateFcns called 

% Hint: edit controls usually have a white background on Windows. 

% See ISPC and COMPUTER. 

if ispc && isequal(get(hObject ,' BackgroundColor ') , 
get (0, ' defaultUicontrolBackgroundColor' )) 

set(hObject,' BackgroundColor' , 'white' ); 

end 

function dzr_Callback(hObject, eventdata, handles) 

% hObject handle to dzr (see GCBO) 

% eventdata reserved - to be defined in a future version of MATLAB 
% handles structure with handles and user data (see GUIDATA) 

% Hints: get (hObject, 'String') returns contents of dzr as text 

% str2double(get(hObject,'String')) returns contents of dzr as a double 

% - Executes during object creation, after setting all properties. 
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function dzr_CreateFcn(hObject, eventdata, handles) 

% hObject handle to dzr (see GCBO) 

% eventdata reserved - to be defined in a future version of MATLAB 
% handles empty - handles not created until after all CreateFcns called 

% Hint: edit controls usually have a white background on Windows. 

% See ISPC and COMPUTER. 

if ispc && isequal(get(hObjectBackgroundColor'), 
get (0, 'defaultUicontrolBackgroundColor' )) 

set(hObject,' BackgroundColor' , 'white' ); 

end 

function dxr_Callback(hObject, eventdata, handles) 

% hObject handle to dxr (see GCBO) 

% eventdata reserved - to be defined in a future version of MATLAB 
% handles structure with handles and user data (see GUIDATA) 

% Hints: get (hObject, 'String') returns contents of dxr as text 

% str2double(get(hObject,'String')) returns contents of dxr as a double 

% - Executes during object creation, after setting all properties. 

function dxr_CreateFcn(hObject, eventdata, handles) 

% hObject handle to dxr (see GCBO) 

% eventdata reserved - to be defined in a future version of MATLAB 
% handles empty - handles not created until after all CreateFcns called 

% Hint: edit controls usually have a white background on Windows. 

% See ISPC and COMPUTER. 

if ispc && isequal(get(hObject,'BackgroundColor'), 
get (0, 'defaultUicontrolBackgroundColor' )) 

set(hObject,' BackgroundColor' , 'white' ); 

end 

function Nzr_Callback(hObject, eventdata, handles) 

% hObject handle to Nzr (see GCBO) 

% eventdata reserved - to be defined in a future version of MATLAB 
% handles structure with handles and user data (see GUIDATA) 

% Hints: get(hObject,'String') returns contents of Nzr as text 

% str2double(get(hObject,'String')) returns contents of Nzr as a double 

% - Executes during object creation, after setting all properties. 

function Nzr_CreateFcn(hObject, eventdata, handles) 

% hObject handle to Nzr (see GCBO) 

% eventdata reserved - to be defined in a future version of MATLAB 
% handles empty - handles not created until after all CreateFcns called 

% Hint: edit controls usually have a white background on Windows. 

% See ISPC and COMPUTER. 

if ispc && isequal(get(hObject,'BackgroundColor'), 
get (0, ' defaultUicontrolBackgroundColor' )) 

set(hObject,' BackgroundColor' , 'white' ); 

end 

function Nxr_Callback(hObject, eventdata, handles) 

% hObject handle to Nxr (see GCBO) 

% eventdata reserved - to be defined in a future version of MATLAB 
% handles structure with handles and user data (see GUIDATA) 

% Hints: get (hObject, 'String') returns contents of Nxr as text 

% str2double(get(hObject,'String')) returns contents of Nxr as a double 
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% - Executes during object creation, after setting all properties. 

function Nxr_CreateFcn(hObject, eventdata, handles) 

% hObject handle to Nxr (see GCBO) 

% eventdata reserved - to be defined in a future version of MATLAB 
% handles empty - handles not created until after all CreateFcns called 

% Hint : edit controls usually have a white background on Windows . 

% See ISPC and COMPUTER. 

if ispc && isequal(get(hObjectBackgroundColor'), 
get ( 0, 'defaultUicontrolBackgroundColor' )) 

set(hObject,' BackgroundColor' , 'white' ); 

end 


% - Executes on slider movement. 

function thetas_slider_Callback(hObject, eventdata, handles) 

% hObject handle to thetas_slider (see GCBO) 

% eventdata reserved - to be defined in a future version of MATLAB 
% handles structure with handles and user data (see GUIDATA) 

% Hints: get (hObject, 'Value') returns position of slider 

% get(hObject,'Min') and get(hObject,'Max') to determine range of slider 

% - Executes on selection change in rx_xdist. 

function rx_xdist_Callback(hObject, eventdata, handles) 

% hObject handle to rx_xdist (see GCBO) 

% eventdata reserved - to be defined in a future version of MATLAB 
% handles structure with handles and user data (see GUIDATA) 

% Hints: contents = cellstr (get(hObject, 'String')) returns rx_xdist contents as 
cell array 

% contents{get(hObject,'Value')} returns selected item from rx_xdist 

% - Executes during object creation, after setting all properties. 

function rx_xdist_CreateFcn(hObject, eventdata, handles) 

% hObject handle to rx_xdist (see GCBO) 

% eventdata reserved - to be defined in a future version of MATLAB 
% handles empty - handles not created until after all CreateFcns called 

% Hint: popupmenu controls usually have a white background on Windows. 

% See ISPC and COMPUTER. 

if ispc && isequal(get(hObject,'BackgroundColor'), 
get (0, ' defaultUicontrolBackgroundColor' )) 

set (hObject, ' BackgroundColor' , 'white' ); 

end 


% - Executes on selection change in rx peddbx. 

function rx_peddbx_Callback(hObject, eventdata, handles) 

% hObject handle to rx_peddbx (see GCBO) 

% eventdata reserved - to be defined in a future version of MATLAB 
% handles structure with handles and user data (see GUIDATA) 

% Hints: contents = cellstr(get(hObject ,' String ')) returns rx_peddbx contents 
as cell array 

% contents{get(hObject ,'Value') } returns selected item from rx_peddbx 

% - Executes during object creation, after setting all properties. 

function rx_peddbx_CreateFcn(hObject, eventdata, handles) 

% hObject handle to rx_peddbx (see GCBO) 

% eventdata reserved - to be defined in a future version of MATLAB 
% handles empty - handles not created until after all CreateFcns called 

% Hint: popupmenu controls usually have a white background on Windows. 
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% See ISPC and COMPUTER. 

if ispc && isequal(get(hObject, 'BackgroundColor' ), 
get (0, 'defaultUicontrolBackgroundColor' )) 

set (hObject, 'BackgroundColor' , 1 white' ); 

end 


% - Executes on selection change in rx_nexpx. 

function rx_nexpx_Callback(hObject, eventdata, handles) 

% hObject handle to rx_nexpx (see GCBO) 

% eventdata reserved - to be defined in a future version of MATLAB 
% handles structure with handles and user data (see GUIDATA) 

% Hints: contents = cellstr(get (hObject,'String')) returns rx_nexpx contents as 
cell array 

% contents { get(hObject ,' Value ') } returns selected item from rx_nexpx 

% - Executes during object creation, after setting all properties. 

function rx_nexpx_CreateFcn(hObject, eventdata, handles) 

% hObject handle to rx_nexpx (see GCBO) 

% eventdata reserved - to be defined in a future version of MATLAB 
% handles empty - handles not created until after all CreateFcns called 

% Hint: popupmenu controls usually have a white background on Windows. 

% See ISPC and COMPUTER. 

if ispc && isequal(get(hObject,'BackgroundColor'), 
get (0, ' defaultUicontrolBackgroundColor' )) 

set(hObject,' BackgroundColor' , 'white' ); 

end 


% - Executes on selection change in rx_zdist. 

function rx_zdist_Callback(hObject, eventdata, handles) 

% hObject handle to rx_zdist (see GCBO) 

% eventdata reserved - to be defined in a future version of MATLAB 
% handles structure with handles and user data (see GUIDATA) 

% Hints: contents = cellstr(get(hObject,'String')) returns rx_zdist contents as 
cell array 

% contents{get(hObject,'Value')} returns selected item from rx_zdist 

% - Executes during object creation, after setting all properties. 

function rx_zdist_CreateFcn(hObject, eventdata, handles) 

% hObject handle to rx_zdist (see GCBO) 

% eventdata reserved - to be defined in a future version of MATLAB 
% handles empty - handles not created until after all CreateFcns called 

% Hint: popupmenu controls usually have a white background on Windows. 

% See ISPC and COMPUTER. 

if ispc && isequal(get(hObject ,' BackgroundColor ') , 
get (0, 'defaultUicontrolBackgroundColor' )) 

set (hObject, ' BackgroundColor' , 'white' ); 

end 


% - Executes on selection change in rx_peddbz . 

function rx_peddbz_Callback(hObject, eventdata, handles) 

% hObject handle to rx_peddbz (see GCBO) 

% eventdata reserved - to be defined in a future version of MATLAB 
% handles structure with handles and user data (see GUIDATA) 

% Hints: contents = cellstr (get(hObject,' String ')) returns rx_peddbz contents 
as cell array 

% contents { get(hObject ,'Value') } returns selected item from rx_peddbz 
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% - Executes during object creation, after setting all properties. 

function rx_peddbz_CreateFcn(hObject, eventdata, handles) 

% hObject handle to rx_peddbz (see GCBO) 

% eventdata reserved - to be defined in a future version of MATLAB 
% handles empty - handles not created until after all CreateFcns called 

% Hint: popupmenu controls usually have a white background on Windows. 

% See ISPC and COMPUTER. 

if ispc && isequal(get(hObjectBackgroundColor'), 
get ( 0, 'defaultUicontrolBackgroundColor' )) 

set(hObject,' BackgroundColor' , 'white' ); 

end 


% - Executes on selection change in rx_nexpz . 

function rx_nexpz_Callback(hObject, eventdata, handles) 

% hObject handle to rx_nexpz (see GCBO) 

% eventdata reserved - to be defined in a future version of MATLAB 
% handles structure with handles and user data (see GUIDATA) 

% Hints: contents = cellstr (get(hObject,'String')) returns rx_nexpz contents as 
cell array 

% contents{get(hObject ,'Value') } returns selected item from rx_nexpz 

% - Executes during object creation, after setting all properties. 

function rx_nexpz_CreateFcn(hObject, eventdata, handles) 

% hObject handle to rx_nexpz (see GCBO) 

% eventdata reserved - to be defined in a future version of MATLAB 
% handles empty - handles not created until after all CreateFcns called 

% Hint: popupmenu controls usually have a white background on Windows. 

% See ISPC and COMPUTER. 

if ispc && isequal(get(hObject ,' BackgroundColor ') , 
get (0, 'defaultUicontrolBackgroundColor' )) 

set (hObject, ' BackgroundColor' , 'white' ); 

end 
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